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during  the  summer  of  1982.  Thanks  are  due  to  Dr.  M.  J.  Werle  for  his  encouragement, 
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ABSTRACT 

.^)  This  work  aims  at  application  in  both  internal  and  external  flows.  An  analytical 
and  computational  study  is  described  which  attempts  ultimately  to  uncover  the  nature  of 
the  flow  structure  surrounding  the  wake  closure  region  of  a  long  recirculatory  eddy 
following  separation  off  an  isolated  aerodynamic  element.  This  problem  area  forms 
a  critical  element  in  the  development  of  a  formally  rational  model  for  a  massive 
separation  eddy  at  large  Reynolds  number.  An  approach  is  developed  here  for  addressing 
the  external  aerodynamic  problem  through  a  family  of  internal,  cascade  flow,  solutions 
with  increasing  gap  spacing.  Each  element  of  this  family  is  formally  correct,  thus 
providing  a  basis  for  systematic  study  of  the  desired  limit  case  of  an  infinite  gap 
spacing. 

The  study  is  in  two  parts.  Part  I  presents  the  structure  for  finite  gap  spacing 
and  the  corresponding  computational  results.  Part  II  discusses  the  possible  flow 
structures  emerging  for  widely  spaced  cascades  and  their  connection  with,  and 
implications  for,  the  external  aerodynamic  flow. 
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PART  I  -  THE  FINITE-SPREAD  CASCADE 


1-1.  INTRODUCTION 


The  study  of  the  fluid  flow  at  high  Reynolds  number  through  a  cascade  of  bluff 
bodies  is  useful,  it  is  believed,  for  two  major  reasons.  First,  the  resultant 
predictions  should  have  some  relevance  to  the  performance  of  the  cascade  itself, 
serving  in  some  sense  as  a  model  for  flow  past  compressor  and  turbine  blades,  in  heat 
exchangers,  and  in  other  interactive  body  configurations.  In  practice,  the  resulting 
breakdown  of  mainly  attached  flow  can  have  severe  consequences  and  so  the  predicted 

drag  and  eddy  lengths,  among  other  things,  are  of  immediate  concern.  The  second 

major  reason  for  considering  the  cascade  configuration  is  that  it  provides  an  alter¬ 
native  means  of  gaining  extra  insight  into  the  properties  of  external  motions. 

The  main  such  properties  are  attained  when  the  cascade  becomes  widely  spaced  apart. 

Since,  as  shown  here,  the  understanding  of  the  separated  flow  structure  produced  by 
the  cascade  configuration  can  be  completed,  this  establishes  a  firm  basis,  therefore, 
for  advancing  gradually  toward  an  understanding  of  the  still  not  fully  known  structure  of 
external  flow  past  a  bluff  body  (whereas  direct  analysis  is  hounded  by  nonuniqueness 
and  nonlinearity).  In  the  latter  case  the  most  puzzling  issue  concerns  reattachment, 
or  the  closure  of  the  recirculating  eddy  or  eddies  behind  the  body;  fortunately,  the 
nature  of  the  reattachment (s)  behind  a  cascade  is  one  of  the  most  important  issues  of 
the  cascade  problem,  too. 

The  present  theoretical  work  supposes  the  Reynolds  number  Re  to  be  large,  the 
fluid  to  be  incompressible  and  its  motion  to  be  steady,  laminar  and  two-dimensional. 

There  is  a  possible  application  [1981,  M.  J.  Werle,  private  communications]  of  the 
theory  to  turbulent  flow  also,  however,  with  regard  to  the  reattachment  process.  The 
whole  cascade  is  taken  to  be  infinite  in  extent,  and  symmetric,  so  that  only  the  flow 
past  the  upper  half  of  one  body  needs  to  be  discussed  provided  the  relevant  lines  of 
symmetry  are  present  (see  Fig.  1).  In  nondimensional  terms  the  chord  of  the  bluff 

body  is  1,  the  given  cascade  spacing  H  is  finite,  of  0(1),  the  uniform  freestream 

c 

speed  is  u  =  1,  and  the  velocities,  corresponding  Cartesian  coordinates,  stream  function, 
flow  speed  and  pressure  are  (u,v),  (x,y),  i}>,  q,  p,  respectively.  The  pressure  level 
far  upstream  is  normalized  to  zero.  The  body  is  "bluff"  in  the  sense  that  its  typical 
thickness  is  allowed  to  be  0(1)  and  its  trailing  edge,  possibly  wedged  or  blunt,  is 
not  cusped. 

The  separated  flow  structure  is  described  in  Section  1-2  below  and  consists  of 

the  body  scale  flow  and  the  long  eddy  or  reattachment  scale  flow,  considered  in  turn. 

The  Kirchhoff  free  streamline  solution  with  a  slowly  moving  eddy  gives  the  dominant 

self-consistent  account  of  the  body  scale  flow,  which  leads  to  a  constant  0(1)  eddy 

width  H  downstream.  The  eddy  is  then  closed,  yielding  reattachment,  on  a  longer  0(Re) 
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length  scale  downstream.  There,  across  the  whole  flow,  the  boundary  layer  equations 
with  unknown  pressure  control  matters,  with  the  viscous  forces  preventing  any  strong 
backward  jet  from  arising,  and  so  overall  self-consistency  is  achieved.  The 
reattachment  problem  depends  only  on  the  ratio  H  =  (see  1)  anc*  its 

numerical  treatment  is  outlined  in  Section  1-3.  The  treatment  involves  a  three- 
regioned  scheme,  designed  to  acknowledge  the  three-layered  form  of  the  starting 
solution  where  the  initial  velocity  profile  is  discontinuous,  and  windward-differencing 
to  accommodate  the  upstream  influence  within  the  reversed  flow.  The  main  results 
obtained  are  presented  for  a  range  of  values  of  H  =  Hw/H  ,  the  nondimensional  eddy 
half  thickness.  Finally,  Section  1-4  presents  our  interim  comments,  preparatory  for 
Part  II,  including  the  connection  with  external  flow  obtained  when  H  tends  to  zero 
and  certain  possibilities  regarding  upstream  influence  in  shear  layers. 

It  is  worth  remarking  here  that  the  present  reattachment  problem  opens  probably 
one  of  the  most  direct  ways  of  approaching  the  puzzle  of  large-scale  external 
reattachment:  see  also  Part  II.  The  way  seems  easier  or  more  direct  than  in  triple¬ 
deck  calculations  for  instance.  At  the  same  time,  the  results  for  finite  cascade 
spacing  could  provide  some  helpful  applications. 
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1-2.  THEORETICAL  FLOW  DEVELOPMENT 


On  the  body  scale,  where  x,  y  are  0(1),  the  separated  flow  structure  is  of  the 
inviscid  Kirchhoff,  free-streamline,  type  to  leading  order  with  appropriate  extensions 
(see  Refs.  1,  2)  to  include  the  effects  of  viscosity.  So  the  following  properties 

(a)  -  (i)  hold. 


(a) 


The  relative  error  in  the  dominant  inviscid  solution  is  0(Re  see 

also  (e)  below.  Thus, 


(u,<f<,p) 


,  v  "1/16  , 

(u  ,i^  ,p  )  +  Re  (u  ,i|»  ,p  )  +  ... 

o  o  o  111 


(1-2.1) 


where  u  ,ii>  ,p  is  the  Kirchhoff  from  (see  sketch  in  Fig.  1). 
o  o  o 

(b)  The  dominant  solution  u  ,  ij>  ,p  then  satisfies  the  condition  of  uniform 

pressure  on  the  separat?ng°streamline  <J>  =  0,  y  =  S  (x)  (Fig.  1),  so  that 

o 


P  =  C  on  y  =  S  (x)  (1-2.2) 

o  o  o 


2 

Here  2C  =  [1-((H  /H,)  /  (H  /H  -1))  ]  is  a  negative  constant,  from  Bernoilli's 

O  C  yy  C  •• 

theorem  and  mass  conservation  considerations,  with  being  defined  in  Fig.  1  and 
in  (g)  below. 


(c)  Within  the  eddy  between  the  separating  streamline  and  the  x-axis,  the 

velocities  are  small,  |q|  <<  1,  and  in  particular  uD  =  iJj  hO,P  =  C. 

It  is  possible  that  in  fact  |q|  is  as  low  as  0(Re“l/2)  in  the  edcly  but 

that  is  the  least  value  possible,  due  to  entrainment,  and  the  principal 

contention  of  the  theory  does  not  yet  need  to  be  made  more  specific  than 

u  =  0 . 
o 

(d)  On  the  body,  between  the  front  stagnation  point  and  the  separation  point, 
the  pressure  gradient  is  favorable,  forcing  automatically  an  attached 
0(Re~l/2)  boundary  layer  there. 


(e)  The  viscous  separation  process  is  local,  of  the  interactive  triple-deck 
form  (see  Refs.  1-3),  and  that  is  responsible  (Ref.  1)  for  the  relatively 
large  correction  required  in  Eq.  (1-2.1). 
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(f)  The  viscous  shear  layer  surrounding  the  separating  streamline  has  thickness 

0(Re  )  end  is  subjected  to  uniform  external  velocities  of  q  %  Hc/ (Hc-Hw) 

and  q  %  0  at  its  upper  and  low  extremes,  respectively.  At  its  start, 
immediately  beyond  the  separation  point,  it  matches  (Ref.  1)  with  the 
interactive  separation  process  of  (e) .  At  its  other  end,  as  x  -*■  °°,  its 
0(Re“l/2)  thickness  grows  like  x^^  as  the  Chapman  similarity  form  is 
acquired. 

(g)  The  eddy  width  SQ(x)  tends  to  a  constant  %  downstream  as  x  -*■  c.f.  the 

unbounded  growth  of  external  flow  (Refs.  1,  2).  Here  Hy  depends  only  on 
the  cascade,  spacing  Hc,  for  a  given  body  shape,  and  is  related  to  the 
principal  drag  Cp  on  the  body. 

(h)  If  Hc  is  large  for  a  given  body  shape,  then  the  dominant  inviscid  solution 
effectively  subdivides.  First,  for  0(1)  values  of  x,y  the  Kirchhoff 

form  for  external  flow  is  retrieved  as  the  presence  of  all  the  other 
bodies  becomes  negligible,  to  leading  order.  So  the  unbounded  growth 
SQ(x)  'v  (4  Cjj/ir)  1/2  x  retrieved  downstream  as  x  +  “  on  that  inner 

subscale.  Second,  however,  the  presence  of  the  0(Hc)  distant  confinement 
reasserts  itself  on  an  outer  subscale  when  x  is  0(HC),  with  SD(x)  then 
being  0(Hc^/2)  for  consistency.  The  adjustment  there  is  similar  to  that 
described  in  the  appendix  of  Ref.  1,  the  earlier  unbounded  growth  is 
gradually  halted  and  as  a  result  the  constant  eddy  width  Hy  is  achieved 
when  x  »  Hc«  See  also  Part  11.  The  property  found  is  that 

1/2 

H  «  H  when  H  is  large.  (1-2.3) 

W  c  c 

A  like  feature  holds  in  effect  if  the  cascade  spacing  Hc  is  kept  fixed  and  the 
body  size  is  reduced  appreciably  instead. 

(i)  Some  inviscid  aspects  of  the  above  structure  are  considered  in  some 
textbooks  (see,  for  example,  Ref.  9). 

The  structure  of  the  body  scale  flow  above  is  believed  to  be  entirely  self- 
consistent  with  regard  to  the  Navier-Stokes  equations.  We  are  left  now  with  the 
task  of  describing  the  ultimate  wake  closure  process  further  downstream,  beyond  the 
body  scale,  as  the  recirculatory  eddy  has  to  close  and  the  motion  then  returns  to 
its  uniform  state  eventually.  The  closure  or  reattachment  takes  place  on  the 
longer  0(Re)  length  scale  in  x  downstream,  since  then  the  shear  layer  thickness  of 
(f)  above  grows  to  0(1).  So  the  expressions 


(u,<|/,p)  =  (0*,**,p*)  +  ...  ,  with  x  =  Re  X*.  y  =  0(1) 


(1-2.4) 


hold  in  the  closure  phase,  with  the  Navier-Stokes  equations  thereby  reducing  to 
the  classical  boundary  layer  equations  (as  in  Re's.  4-6),  for  X  >  0, 
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U 


* 


3tJj* 

3y 


3X* 


3^*  3U* 


* 

3X 


3y 


dP*  32U* 
dX*  3y2 


(1-2. 5a) 


with  P  dependent  only  on  X  .  The  formal  relative  error  involved  in  Eq.  (1-2. 5a) 
is  remarkably  small,  of  the  order  Re-2.  The  starting  conditions  required  are  those 
of  a  step-like  velocity  profile  initially: 


"k  k 

at  X  =  0+,  U 


0  for  0  <  y  <  H 
-  3  W 

V'W  f°r  "w < y  i  Hc> 


<1-2. 5b) 


to  match  with  the  slow  eddy  velocity  and  the  Kirchhoff  solution  closer  to  the  body. 
The  boundary  conditions  in  y  are 


ij;*  =  3U*/3y  =  0  at  y  =  0 


(1-2. 5c) 


4i*  =  H  ,  3U*/3y  =  0  at  y  =  H  , 
c  c 


(1-2. 5d) 


k  k  k 

for  X  >  0,  from  the  cascade  symmetry.  The  induced  pressure  P  (X  )  is  unknown. 

The  basic  closure  problem  is  totally  contained  in  Eqs.  (I-2.5a-d),  which 
require  a  r,->merical  solution.  The  present  standpoint  is  that  the  overall  theory 
is  complete  since  the  presence  of  significant  viscous  forces,  in  Eq.  (1-2. 5a), 
during  the  closure  is  always  sufficient  to  allow  the  constraint  of  Eq.  (1-2. 5b) 
to  be  met.  This  constraint  is  necessary  to  prevent  a  strong  backward  jet  from 
entering  the  body  scale  flow  and  provided  that  is  satisfied  the  earlier  structure 
then  remains  intact.  The  crucial  action  of  viscosity  here  is  due  to  the  flow 
confinement  caused  by  the  cascade  symmetry,  which  in  turn  restricts  the  lateral 
distance  y  to  0(1)  values.  This  is  in  direct  contrast  to  the  expensive  dimensions 
reached  by  external  eddies  (according  to  the  same  theory)  where  the  influence  of 
viscosity  is  correspondingly  much  less  clear  and  where  in  consequence  our  under¬ 
standing  has  still  to  be  completed  (see  also  Section  1-4  and  Part  II  for  further 
discussion  of  the  point).  Five  main  features  of  the  current  closure  problem  can 
be  delineated,  in  ( 1 )—  (5 )  below,  in  terms  of  the  renormalized  form  of  the 
governing  equations 
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0  , 
3Y  ’ 

3U  3j)  3U  dP*  3U2 

U  3X  "  3X  3Y  "  "  dX  3y2 

(1-2 . 6a) 

(  0  for  0  <  Y  <  H 

(1-2. 6b) 

at  X  =  0  +, 

P  =  C  ,  U  =  J 

( 1/ ( 1-H)  for  H  <  Y  <  1 

(1-2. 6c) 

ip  = 

3U/3Y  =  0  at  Y  =  0 

(1-2. 6d) 

-  l. 

3U/3Y  =  0  at  Y  =  1 

(1-2. 6e) 

obtained  by  setting,  in  Eq.  (1-2. 5ad), 


(U*,**,P*,X*,y)  =  (U,  Hc^,  P*,  Hc2X,  H  Y), 


(1-2.7) 


with 


“  5  V“c  • 


(1-2.8) 


defining  the  governing  parameter  H. 

(1)  The  solution  depends  only  on  the  parameter  H,  the  ratio  of  the  incoming 
eddy  width  to  the  cascade  spacing. 

(2)  Because  of  the  initial  discontinuous  profile  in  Eqs .  (1-2. 6b, c)  the 
solution  is  singular  as  X  ->  0  +,  with  the  Chapman  similarity  form  being 
recovered  then.  Thus,  near  Y  =  H,  for  small  X, 

1/2  1/2 
ip  X  f  (n)  +  ...  ,  with  n  =  (Y-H)/X  ,  (1-2. 9a) 

c 

where  from  Eqs.  (I-2.6a-c)  f  satisfies  the  Blasius  equation  but  with 
shear  layer  conditions: 

f  +  ;f  f  ”  =  0,  f  ’(»)  =  1/(1-H),  f  '(-«)  =  0  (1-2. 9b) 

c  2  c  c  c  c 


The  solution  of  Eq.  (1-2. 9b)  gives 


UTRC83-13 


f  (-») 
c 


-<  (1-H) 


-1/2 


(1-2. 9c) 


where  <  =  1.24,  although  there  is  an  uncertainty  in  the  origin  for 
r\  (see  also  point  (3)  below) . 

(3)  The  initial  development  of  the  solution  for  small  X  is  three-layered  but 

is  not  uniquely  determined  by  local  considerations.  Generally,  Eq .  (1-2. 9b) 
yields  the  property 


f  a,  (n-q  )/  (l-H)  as  n  -*•  00  (l-2.9d) 

c  o 


where  n  is  an  unknown  constant  origin  shift.  As  a  result,  above  the 
Chapman  layer 


/  Y"H  \  1/2 

*  "  \  i_H  /  +  X  *1/2  (Y)  +  *•' 

for  H  <  Y  <1,  with  Eqs.  (1-2. 6a, e)  requiring 

*1/2  (Y)  ’  -(1-H)  Pl/2  (1-Y)i  P*  '  Co  '  Pl/2  1(1/2  + 


(1-2. 10a) 


(1-2. 10b) 


Hence,  from  merging  with  Eq.  (1-2. 9d),  we  have 


P  =  +  n  /(l-H)3  . 
1/2  o 


(1-2. 10c) 


Then,  from  the  form  of  P  (X)  in  Eq.  (1-2. 10b),  the  solution  below  the 
Chapman  layer  has  the  expansion 

=  Xl/4  G(Y)  +  ...  (1-2. 11a) 

for  0  <  Y  <  H,  since  Eq.  (1-2. 6b)  holds.  Substituting  Eqs.  (1-2. 11a)  and 
(1-2. 10b)  into  Eqs.  (1-2. 6a, d)  and  imposing  G(H-)  =  0  in  view  of' 

Eq.  (1-2. 9a, c),  it  is  found  that 

G(Y)  =  +  (2  P,.„)1/2  —  sin  (nir  Y/H)  (1-2. lib) 

“  1/2  nn 
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where  n  is  a  positive  integer.  Physical  sense  and  other  considerations 
(see  Ref.  5)  suggest  that  n  =  1  here  and  then  the  -  sign  is  the  appropriate 
one.  The  constant  n0  remains  undetermined  throughout;  moreover,  an  even 
larger  origin  shift  in  Y  itself  is  possible  in  Eq.(I-2.9a).  If  n0  “  0 
then  nonuniqueness  can  still  arise  in  higher  order  terms  of  the  three¬ 
layered  expansion  Eqs.  (1-2. 9a),  (1-2. 10a)  and  (1-2. 11a).  A  special 
case  however  is  that  with  n  =0  and  with  the  flow  below  the  Chapman 
layer  being  forced  by  the  entrainment  conditions,  Eq.  (1-2. 9c),  giving 


ip  =  — k  X1/2  (1-H)~1/2  Y/H  + 


there,  which  requires  the  pressure  expansion  to  be 


*  Ik 

p  =  c  -  - - rx+  ...  . 

(l-H)H 


Along  with  Eq.  (1-2. 12a, b)  we  then  have 


(1-2. 12a) 


(1-2. 12b) 


■  (s-) 


* ... 

2  H2 


(1-2. 12c) 


above  the  Chapman  layer. 

It  follows  from  point  (3)  above  that  reversed  flow  is  present  immediately 
in  X  >  0,  below  the  Chapman  layer.  This  feature  also  seems  in  line  with 
the  nonuniqueness  arising  above  since  the  governing  equations  are  locally 
parabolic  in  the  flow  direction,  yielding  upstream  influence  throughout 
0  <  X  <  Xreatt  where  X  =  Xreatt  is  the  closure  position. 


(5)  Downstream,  as  X  +  the  uniform  state  is  recovered  with  U 
for  all  X  the  integrated  momentum  balance 

1 

P*(X)  +  J  U2  (X, Y)  dY  =  Cq  +  1/  (1-H) 


holds,  from  Eq.  (I-2.6a-3).  Therefore  the  asymptotes 


1 .  But 


(1-2.13) 
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1  /  H  \2 

U  -*■  1,  P  -*• 

“  2  (l 15/  3SX  — 

(1-2.14) 

give  the  far  downstream 

behavior. 

Henceforth  we  deal  with 

the  modified  pressure  P 

£ 

P  -  C  ,  tor  convenience 
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1-3.  NUMERICAL  TREATMENT  OF  CLOSURE/REATTACHMENT 


The  aspects  of  points  (2),  (3)  and  (4)  noted  above  -  the  singularity  due  to 
the  discontinuous  initial  profile,  the  initial  three-layered  structure  and  the 
upstream  influence  in  the  reversed  flow  -  would  seem  to  be  so  fundamental  to  the 
closure  problem  that  a  satisfactory  numerical  treatment  must  acknowledge  them 
ultimately.  Perhaps  the  most  important  aspect  however  is  that  of  point  (4), 
which  requires  a  numerical  treatment  to  be  not  solely  forward -marching .  To 
accommodate  reversed  flow  we  used  a  windward-differencing  approach  (see  e.g.. 

Ref.  7).  This  was  used  within  three  progressively  more  refined  schemes  I,  II, 

III:  scheme  I  adopted  uniform  differencing  in  X,  Y;  scheme  II  adopted  differencing 

in  coordinates  sketched  in  the  X,  Y  directions;  and  scheme  III  applied  a  three- 
region  technique  (Refs.  5  and  8)  to  allow  also  for  the  features  of  points  (2), 

(3)  above.  All  three  schemes  proved  convergent  with  the  windward  differences 
included,  but  we  focus  henceforth  on  the  nominally  most  accurate  one,  scheme  III. 

The  computational  domain  was  split  into  three  regions  A-C,  in  line  with 
point  (3)  above.  Region  A  was  taken  to  be  in  effect  below  the  middle  Chapman 
region  B  and  Region  C  was  taken  above  it  (see  Fig.  2).  A  new  variable  £,  =  X^/2 
is  defined  in  A,  with  ip  =  g(5,Y),  u  =  c(5,Y),  3U/3Y  =  e(5,Y),  leaving  the  first 
order  governing  equations  as 


c=-^*  e  =  —  c  —  _-5&e  =  —  —  +2E  — 

3Y  ’  3Y*  3?  H  d?  4  3Y 


(1-3. la) 


from  Eq.  (1-2. 6a),  with  the  boundary  conditions  from  Eq.  (1-2. 6d)  as 


g=e=0atY=0  (1-3. lb) 

In  B,  the  variables  were  written  tp  =  5f(5,n),  U  =  s(5,n),  3U/3Y  =  5  t(5.n) 

because  of  Eq.  (1-2. 9a),  so  that  here  the  equations  become 


r 


t  T7  +  dP/d5 


(1-3.2) 


Then  in  C  the  solution  has  the  same  structure  as  in  A,  giving  Eq .  (1-3. la)  again 
but  with  the  boundary  conditions  from  Eq .  (1-2. 6e)  as 


g  =  1,  e  =  0  at  Y  =  1. 


(1-3.3) 


: 


UTRC83-13 


At  any  station  £  the  three  regions  are  inter-linked  (Ref.  8)  by  ensuring  that 
4>,  U,  3U/3Y  are  all  continuous  at  the  upper  and  lower  boundaries  n  =  n±<*>,  Y  =  H 
+  £,  of  region  B.  Uniform  differencing  in  n  or  Y  as  appropriate  is  applied 

to  Eqs.  (1-3. la)  and  (1-3.2)  and,  with  the  4  conditions  of  Eqs.  (1-3. lb)  and  (1-3.3) 
and  the  6  inter-linking  relations,  sufficient  information  exists  to  determine  the 
solution  at  that  £  station  upon  applying  (Ref.  8)  Newton  iteration,  matrix  inver¬ 
sion  and  knowing  the  solution  at  the  stations  £  ±  A£.  The  differencing  in  the  £ 
direction  is  allowed  to  account  for  the  local  flow  direction  however,  consistent  with 
the  reversed  local  influence.  If  U  >  0  then  the  value  of  U  is  made  to  depend  on  the 
value  at  £-A£,  as  obtained  with  differences  for  U3U/3X  centered  at  £-A£/2.  If  U  <  0 
then  the  latest  known  values  at  £,  £  +  A£  are  used  instead  with  a  forward  difference 
for  U3U/3X.  The  starting  solution  at  £  =  0  is  set  in  the  obvious  way  for  regions  A, 

C  whereas  that  for  region  B  requires  solving  Eq.  (1-3.2)  at  £  =  0,  i.e.,  solving  the 
Chapman  form  of  Eq.  (1-2. 9b)  there. 

The  domain  is  swept  and  updated  N  times,  each  sweep  passing  from  £  =  0  to 
a  sufficiently  far  downstream  location  £  =  £<*>  where  Eq.  (1-2.14)  is  effectively 
obtained,  until  the  system  converges.  The  convergence  rate  was  accelerated 
in  later  trials  by  sweeping  back  and  forth,  marching  from  £  =  0  to  £  =  £oo, 
then  marching  back  from  £«  to  0,  and  so  on,  which  admits  the  upstream  influence 
faster.  The  upstream  influence  enters  the  solution  sweep-by-sweep  through  the 
windward-differencing  of  the  term  U3U/3X;  note  that  the  term  in  3i^/ 3X  is  not 
windward-differenced,  thus  allowing  Eq.  (1-2. 6c)  to  have  its  necessary  influence 
also.  Various  appropriate  values  of  njrD0,A£,  of  the  steplengths  An,  AY  in  n,  Y 
and  of  the  Newton  iterative  tolerance  were  taken  and  tests  on  the  accuracy  of  the 
results  were  made.  Typically,  satisfactory  overall  convergence  is  achieved  when 
N  *  10  for  moderate  values  of  H,  although  this  increases  for  smaller  values  of  H. 
Further  comments  are  given  in  the  next  section. 
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1-4.  RESULTS  AND  DISCUSSION 


The  main  results  obtained  for  various  values  of  H  between  0.9  and  0.2 
are  shown  in  Fig.  3-5.  In  most  cases  the  constant  nQ  of  Eq .  (1-2. 9d)  was  chosen 
to  be  zero,  with  nonzero  values  having  but  a  very  small  overall  effect. 

Figure  3(a)  provides  an  overview  of  the  flow  behavior  by  monitoring  the 
center  line  velocity,  and  the  top  edge  velocity,  U-j.,  for  relatively  large 

values  of  H,  corresponding  to  small  cascade  gap  spacing.  In  all  cases  considered 
the  closure  process  is  seen  to  occur  quite  rapidly,  with  the  normalized  position 
of  the  centerline  flow  reversal  always  within  half  a  gap  width  of  the  initial 
plane. 

It  is  immediately  striking  from  Figs.  3(a)  and  3(b)  that  as  H  decreases 

the  closure  distance  X  shrinks  very  fast  and  that  the  minimum  centerline 

reatt 

velocity  Um^n  decreases  but  less  rapidly,  while  the  flow  properties  near  closure 
become  more  abrupt.  For  larger  values  of  H  the  features  are  rather  milder,  in 
keeping  with  an  asymptotic  analysis  which  shows  that 


reatt 


(1-H) 


-1/2 


U  .  “  (1-H)  1/2.  asH+1  -. 

mm  * 


(1-4.1) 


The  extreme  smallness  of  Xreatt  once  H  is  less  than  about  1/3  is  more  clearly 
shown  in  Fig.  3(b)  in  terms  of  the  centerline  velocity  behavior.  This  aspect  of 
the  flow  adds  to  the  computational  task  of  course  and  in  practice  it  was  found  to 
limit  the  applicability  of  the  first  two  numerical  schemes  I,  II  noted  in  Section 
3.  It  also  emphasized  the  overriding  need  for  accuracy  in  the  calculations,  for 
checks  to  be  made  on  the  effects  of  mesh  size  (see  Figs.)  and  for  the  proper 
notice  to  be  taken  of  the  significant  upstream  influence.  This  last  is  more 
clearly  demonstrated  in  Fig.  3(c)  which  compares  the  centerline  velocity  distri¬ 
butions  with  and  without*  upwind  differencing  in  the  reversed  flow  region. 

The  computational  method  and  program  appeared  to  work  satisfactorily  for 
moderate  values  of  H  but  difficulties  became  pronounced  as  H  decreased.  This  was 
due  mainly  to  the  relatively  fast  change  in  the  solution  near  wake  closure  then 
and  the  resultant  switching  between  forward  and  backward  differences  there,  which 
led  to  characteristic  oscillations  from  sweep-to-sweep  of  the  flowfield.  Incorpo¬ 
rating  relaxation  in  the  updating  per  sweep,  or  fixing  the  choice  of  directional 
differencing  for  several  sweeps,  or  refining/stretching  the  mesh  near  wake  closure, 
all  alleviated  the  problem  somewhat  but  not  permanently,  and  the  same  was  found  for 
other  modifications  made  in  view  of  certain  analytical  features  for  small  H  suggested 


*This  result  is  actually  achieved  in  the  first  sweep  of  the  iterative  algorithm  by 
zeroing  out  the  longitudinal  convective  term  for  reversed  flow,  an  idea  originally 
Introduced  by  Reyhner  and  Flugge-Lotz  (Ref.  16). 
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by  Part  II  below.  The  flow  properties  are  better  viewed  perhaps  in  the  streamline 
patterns  and  velocity  profiles  of  Figs.  4(a)~(e),  showing  the  increasingly  rapid 
turning  of  the  flow  toward  the  wake  closure  point  as  H  decreases  from  0.9  to 
0.25.  Obviously  as  H  decreases  the  computational  difficulty  is  aggravated,  raising 
serious  concerns  for  the  proper  numerical  modelling  of  the  convective  processes 
then. 


The  pressure  response  P  and  closure  distance  Xreatt  are  given  in  Fig.  5  for 
various  values  of  H.  It  is  noteworthy  that  P(Xreatt)  is  close  to  H  and  Xreatt 
is  closely  proportional  to  [see  log-log  plot]  at  the  smaller  values  of  H:  c.f. 
Part  II. 

As  regards  the  hope  that  an  inkling  of  the  unknown  structure  for  external  flow 
past  a  bluff  body  should  emerge  as  H  becomes  small,  the  major  questions  perhaps 
concern  the  bahavior  of  Xreatt  and  Xmin  (the  position  at  which  U  =  Umin) ,  as 
H  ->  0.  We  might  expect  both  to  tend  to  zero  rapidly  then,  in  line  with  the  pre¬ 
sent  numerical  results  (Figs.  3-5),  but  what  matters  is  the  rate  of  these  trends*. 
Suppose,  for  instance,  that  Xreatt  1  as  H  ->  0.  Then  the  unsealed  closure 
distance,  which  is  x  =  Re  Hc^  Xreatt  from  Eqs.  (1-2.4)  and  (1-2.7),  remains 

0(Re)  as  Hc  +  ®  with  tfae  body  fixed,  since  Eqs.  (1-2.8)  and  (1-2.3)  link  H  to  Hc 
then.  Accordingly,  the  eddy  length  remains  0(Re)  and  that  leads  directly  to  the 
extended  Kirchhoff  account  of  external  bluff  body  flow.  If  instead  Xreatt  <* 
as  H  -*■  0,  however,  then  xreatt  a  Re  H  ^  «  Re  Hc~  for  H,  >>  1  by  similar 

reasoning.  So  then  Xreatt  falls  to  0(S)  as  Hc  is  increased  to  0(Re)  (the  order 
required  for  external  flow  properties  to  be  most  relevant:  see  Ref.  1)  and  that 
points  strongly  to  the  alternative  Prandtl-Batchelor  account  for  external  flows. 

See  also  Part  II. 

Again,  therefore,  there  is  a  call  for  sufficient  accuracy  in  any  further 
calculations  at  even  smaller  values  of  H,  to  help  settle  the  matter  of  the 
external  closure  properties. 

In  the  meantime,  it  is  interesting  that,  on  the  body  scale,  the  Kirchhoff 
solution  for  external  flow  does  emerge  as  Hc  becomes  large  and  the  maximum  eddy 
width  increases  without  bound  then.  It  would  be  remarkable  if  this  process 
could  be  reversed,  at  greater  cascade  spacing,  by  the  change  in  closure  properties 
considered  above.  This  issue  is  amplified  in  Part  II. 

Various  features  of  the  limit  structure  for  the  closure  problem  as  H  +  0 
are  considered  in  Part  II.  A  relatively  thin  detached  shear  layer  tends  to  be 
promoted  from  the  Chapman  region  then  but  the  means  of  reversal  of  its  entrained 
fluid  (where  \p  <  0)  during  closure  at  X  =  Xreatt  has  to  be  established,  given  the 
requirement  of  effectively  zero  backflow  at  X  =  0+.  The  means  may  involve  a 


A  less  sensitive  quantity  to  examine  may  be  the  ratio  Xm^n/Xreatt  which  defines 
to  some  extent  how  local  the  reattachment  process  is  for  a  given  H.  However,  any 
indication  of  the  correct  limiting  form  for  H  ->  0  is  useful. 


15 


UTRC83-13 


localized  strong  interaction  near  the  closure  point,  in  view  of  the  above  discussion, 
which  requires  the  detached  shear  layer  to  admit  some  upstream  influence  ahead  of 
the  closure.  One  kind  of  upstream  influence  studied  is  discussed  in  detail  in 
Appendix  A.  There  it  is  shown  that  upstream  influence  ij3  possible  (see  Eqs. 

A-4a-f)  and  an  interaction  between  the  oncoming  shear  layer  and  the  motion  beneath 
it  can  start.  The  interaction  is  found  to  become  nonlinear  as  X  becomes  finite 
and  it  involves  the  setting  up  (due  to  Ug)  of  a  significant  motion  beneath  the 
shear  layer.  The  ultimate  properties  of  this  interaction  could  be  significant  and 
might  enable  the  shear  layer  to  rid  itself  of  all  the  entrained  fluid  eventually, 
in  readiness  for  the  closure,  for  instance.  Again  further  study  of  the  consequences 
of  the  upstream  interaction  of  Eqs.  (A-l  -  A-4g)  below  should  be  well  worthwhile. 
Appendices  B  and  C  summarize  other  means  of  upstream  influence.  Like  Appendix  A 
they  bring  viscosity  into  play,  the  importance  of  which  is  stressed  in  Part  II. 

They  also  have  certain  repercussions  for  separated  flow  stability.  Another  form 
of  upstream  influence  which  is  currently  under  consideration  and  may  be  relevant 
is  that  occurring  in  separated  flow  as  studied  in  Ref.  10,  on  the  airfoil  trailing 
edge  stall/separation  problem. 

Whatever  the  limiting  form  of  the  present  closure  problem  is  as  H  -+•  0,  all 
the  above  serves  to  emphasize  a  major  point:  the  resolution  of  the  puzzle  concerning 
large-scale  external  closure  is  contained  within  the  standard  boundary  layer  equa¬ 
tions.  We  believe  that  still  further  examination  on  both  the  numerical  and  the 
analytical  sides  of  the  present  basic  topic  could  be  very  rewarding,  therefore. 

Part  II  below  addresses  the  limit  problem  for  H  0  and  its  connection  with 
external  flow. 
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PART  II:  THE  WIDE-SPREAD  CASCADE  AND  EXTERNAL  FLOW 


II-l.  INTRODUCTION 

Given  the  basic  structure,  results  and  guidelines  for  cascades  of  finite 
spread  H-^  in  Part  I,  our  concern  henceforth  is  with  understanding  the  wake  closure 
in  a  wide-spread  cascade  (H-*0)  and  then  its  connection  with  the  external  motion 
past  an  isolated  bluff  body.  There  are  two  main  contenders  for  the  wake  flow  des¬ 
cription,  the  Kirchhoff  free-streamline  solution  (extended)  and  the  Prandtl-Batchelor 
model.  The  latter  tends  to  flounder  for  any  slender  wake,  however,  due  to  the  lack 
of  a  sustained  pressure  gradient,  whereas  Kirchhoff 's  description  already  works 
apparently  for  the  long  slender  wake  of  the  finite-spread  cascade.  So  the  essential 
matter  is  whether  or  not  the  Kirchhoff  solution  continues  to  hold  good  for  wide¬ 
spread  cascades  also. 

Wide-spread  cascades  require  analysis  of  just  the  boundary  layer  equations, 
from  Part  I,  for  wake  closure  as  H-K),  like  the  flow  beyond  a  step-like  trailing 
edge.  This  is  started  in  Section  II-2  below,  which  considers  the  underlying  scales 
and  sublayers  implied.  The  Prandtl-Batchelor  model  can  soon  be  ruled  out.  Indeed, 
near  the  wake  centerline  Y  =  0  all  inviscid  mechanics  fails  to  maintain  consistency 
at  wake  closure,  thus  calling  for  significant  viscous  action.  Section  II-2  also 
stresses  the  0(H)  pressure  rise  P  and  0(H  )  length  scale  in  X  suggested  by  the 
Part  I  calculations  although  the  need  for  wake  closure  then  is  not  definite.  On 
the  contrary.  Section  II-3  shows  that  self-consistency  is  apparently  maintained 
by  wake  closure  on  the  0(1)  scale  in  X,  where  the  required  viscous  action  comes  into 
play  accompanied  by  a  significant  pressure  gradient.  Alternative  accounts  present 
themselves,  however,  due  to  the  nonlinear,  recirculating,  nature  of  the  earlier 
eddy  motion,  as  discussed  in  Section  II-4.  The  Section  II-4  findings  lead  on  to  the 
work  in  Sections  II-5  and  II-8  in  fact. 

Section  II-5  considers  viscous  wake  closure  occurring  with  negligible  pressure 
gradient,  unlike  the  closure  in  Section  II-3.  Similarity  wake  solutions  are  found 
which  yield  algebraic  dependence  at  the  edge  of  a  shear  layer  on  the  verge  of  closure. 
Other  similarity  solutions  exist,  one  showing  that  Goldstein’s  (Ref.  11)  form  of 
the  near-wake  of  an  aligned  flat  plate  is  not  unique:  see  also  Ref.  17. 

Another  spin-off  from  Section  II-5  concerns  reattachment  onto  a  solid  surface. 

The  reattachment  may  be  achievable  with  negligible  pressure  response,  a  desirable 
property,  and  secondary  separation  can  be  induced  further  upstream.  Applications 
of  wide  interest  exist,  to  flow  past  ramps,  over  bluff  bodies  with  splitter  plates, 
or  in  wind-tunnels,  and  these  matters  are  considered  in  Section  II-6. 

In  Section  II-7,  external  flow  and  its  relation  to  the  wide-spread  cascade  are 
studied.  The  linear  dependence  of  eddy  length  on  Reynolds  number  Re,  as  indicated 
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by  previous  studies  (Refs.  4-6),  is  supported  by  an  examination  of  certain 
exponentially  small  terms  provoked  by  Kirchhoff's  solution.  In  turn,  the  earlier 
limit  forms  for  H-+0  all  support  the  emergence  of  Kirchhoff's  solution  for  external 
flows,  on  the  body  scale. 

A  proposed  mechanism  for  the  subsequent  lengthy  wake  closure  is  provided  by 
non-entraining  viscous  shear  layers,  implied  by  Sections  II-4,  II-5  and  addressed 
in  Section  II-8.  Their  four  main  attributes  are:  first,  they  yield  algebraic 
decay  of  the  velocity  in  the  lower  reaches  of  the  shear  layer;  second,  they  show 
that  Chapman's  solution  for  a  shear  layer  far  downstream  is  nonunique;  third,  they 
avoid  the  severe  difficulties  associated  with  supplying  entrainment  to  classical 
shear  layers;  fourth,  they  suggest  a  means  for  massive  eddy  closure  to  take  place 
without  provoking  a  strong  backward  jet  along  the  wake  centerline.  Final  comments 
are  presented  in  Section  II-9. 
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I I— 2 .  ON  LENGTH  SCALES  FOR  THE  WIDE-SPREAD  CASCADE 


From  Part  I,  the  wake  of  the  cascade  with  finite  spread  is  governed  by  the 
boundary  layer  equations 


subject  to 


2 

3tl>  3U  3U  ,  ,  3  U 

U  “  3Y  ’  U  3X  “  3X  3Y  "  ~  P  X)  +  ^ 2 


(II-2 . la) 


3U 

f=-  =  0atY  =  0,X>0 


( X  I  — 2 . lb) 


3U 

=°atY=1,  X>0 


(II-2.1c) 


U  = 


i(l-H)  1  for  H  <  Y  -  1 


0  for  0  -  Y  <  H 


at  X  =  0. 


( I I— 2 . Id) 


Our  present  concern  is  with  the  flow  behavior  for  wide-spread  cascades,  where  H-*0. 

The  following  tentative  arguments  (see  also  Fig.  6)  may  be  made  then. 

At  sufficiently  small  distances  X  ^  A  from  the  start  a  thin  viscous  shear  layer 
is  present,  originating  from  the  Chapman  form  and  positioned  a  distance  0(H)  (at 
Y  =  HS(X/A)  +. . . ,  say)  from  the  center  line  Y  =  0,  in  view  of  (II-2.1d).  The  charac¬ 
teristic  velocity,  mass  flux  and  width  of  the  shear  layer  are  of  the  orders  1,  A^^, 
A1  in  turn.  The  pressure  force  has  negligible  effect  here  since,  above  the  shear 
layer,  on  a  Y  scale  of  0(1),  inviscid  linearized  properties  hold,  yielding  the 
pressure-displacement  law  of  one-dimensional  channel  flow,  P  =  H  -  HS,  for  mass 
conservation. 

Hence,  on  this  scale  we  have 


m 


m 


X  =  AX,  P  =  HP (X)  +  .  . . 


( 1 1  —  2 . 2a) 


m 


-  -  1/2 

with  P(0)  =  0  and  in  the  shear  layer  I,  where  Y  =  HS(X)  +  A  Yq,  (<p,  U)  = 
(/\l/2  '  t  the  pressure-free  boundary  layer  equations 


U 


o 


3<1  3U0 

o  _  _ 

3Y  ’  U°  3X 
o 


2 

'lip 

3  U 

o 

0 

o 

•X 

3Y 

„  2 

o 

3Y 

o 


( 1 1-2 . 2b) 


<9 

“•  .1 
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apply.  Here  and  below  the  Y-coordinate  could  be  displaced  further  If  necessary 
because  of  Prandtl’s  transformation.  The  constraints  here  are  the  familiar  ones 
for  detached  shear  layers, 


U  (X.->  -  1,  U  (X.-)  -  0,  0  (0,Y  ) 

o  o  o  o 


1  for  Y  >  0 
o 

0  for  Y  <  0 
o 


(II-2.2c) 


essentially  from  the  starting  form  (II-2.1d).  The  Chapman  singularity 


-1/2 
\l>  -v-  X 

o 


fc(n),  X-»0,  n 


y  nin 

o 


(11-2. 2d) 


holds  at  the  start.  Above  the  shear  layer  where  Y  is  0(1)  (region  II)  we  have 


(♦,U)  -  [Y  +  H(l-P) (Y-l) ,  1  +  H(l-P) ]  +  0(H)  (II-2.3) 

for  mass_and  momentum  conservation,  from  (II-2.1a,  c,  d) .  Matching  with  layer  I  as 
Y  -*•  H  S(X)  +  then  gives 


P(X)  -  1  -  S(X) , 


(11-2.4) 


as  anticipated  above.  Further  terms  in  the  expansions  of  regions  I,  II,  and  III 
below,  may  be  written  down  in  principle  but  the  balancing  of  the  leading  order  con¬ 
tributions  sets  the  major  task. 

A  complete  solution  for  the  shear  layer  of  (II-2.2b-d)  is  the  Chapman  form 
(II-2.2d)  but  holding  for  all  X>0.  We  take  that  as  the  solution  for  now  (although 
it  is  found  later  to  be  nonunique:  see  Section  II-8)and  follow  its  implications. 

The  main  point  then  is  that,  since  fc(-°°)  =  -  tc  <  0,  a  finite  fraction  of  the 
viscous  shear  layer  I  contains  fluid  entrained  from  below  I,  with  <  0.  This 
fluid  must  be  ejected  further  downstream  before  reattachment/closure  of  the  wake 
eddy  or  eddies  can  occur.  The  viscous  entrainment  forces  the  condition 


* 


-1/2  1/2 
k  X  A 


(-*E)  as  Y  HS(X)  -  , 


(II-2.5) 


which  causes  the  fluid  below  I  to  be  set  into  motion. 
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3 

Immediately  the  length  scale  A  =  0  (H  )  suggests  itself.  For,  as  P  is 
typically  0(H)  from  (II-2.2a),  then  predominantly  inviscid  mechanics  in  the  0(H) 
thick  layer  III  (occupying  0  <  Y  <  HS(X))  beneath  I  implies  that  U  is  0(H^'2)  there, 

O  /  O 

corresponding  to  an  iJj  of  0(H  '  ),  due  to  the  pressure  response.  So  a  comparison 
with  the  viscous  condition  (11-2*5) indicates  that  A  =  0(HJ)  may  represent  a 
crucial  stage  in  the  flow  development,  as  indeed  the  calculations  of  Part  I  suggest 
also.  On  a  shorter  scale  A  <<  H^  the  entrainment  is  negligible  compared  with  a 

*i/o 

pressure- induced  ^  of  0(H  '  )  in  layer  III,  so  that  such  a  flow  in  layer  III,  if 
nontrivial,  must  recirculate  between  I  and  the  center  line.  This  is  not  implausible, 
but  it  leads  to  no  clear  sign  of  eddy  closure  or  anything  of  significance.  For 
Bernouilli’s  law  holding  in  layer  III  does  not  allow  positive  pressure  P,  from  the 
center  line  flow,  therefore  S(X)  >  1  from  (II-2.4)  and  that  excludes  the  closure 
requirement  S  -*■  0.  The  same  objection  must  be  overcome  for  all  inviscid  analysis 
along  the  centerline. 

In  passing  we  note  four  asides.  First,  a  reasoning  similar  to  the  above  tends 
to  rule  out  the  Prandtl-Batchelor  proposal  of  a  strong  mainly  inviscid  recirculating 
eddy  of  uniform  vorticity,  as  follows.  If  such  an  eddy  is  thin,  then  it  induces 
little  pressure  change,  from  (II-2.4),  so  that  the  bounding  forward  shear  layer 
cannot  be  forced  back  at  closure.  On  the  other  hand,  if  the  eddy  is  thick,  with 
Y  of  0(1)  [0  <  Y  <  s(X),  say],  then  between  it  and  Y  =  1  the  inviscid  properties 
holding  allow  no  significant  vorticity  92^/9Y  because  of  the  uniform  starting 
profile  in  (II-2.1d),  and  mass  and  momentum  conservation  then  give  the  nonlinear 
cov  cerpart  of  (II-2.4),  P  =  (s2/2  -  s)/(l  -  s)2  since  s(0)  =  0  to  leading  order; 
however,  constant  vorticity  32’]//9Y2  =  n  inside  the  eddy  yields  ij;  =  QY(Y-s)/2, 

U  =  Q(Y  -  s/2),  P  =  -  Q2  s2/8  on  integration;  so  the  pressure  match  demands  the  final 
result  s(X)  =  0.  Both  cases  lead  to  failure,  therefore.  The  above  requires  the 
length  scale  A  to  be  small,  incidentally,  study  of  the  0(1)  scale  in  X  being  deferred 
to  Section  II-3  below.  Second,  viscous  forces  in  the  lowest  layer  III  remain 
secondary  effects  provided  |U9U/9X|  =  0(HA  ^)  is  larger  than  |9  U/9Y2|  =  0(H-^2), 

i.e.,  for  A  <<  H^  ,  given  the  proposed  scalings  of  (II-2.2a)  -  (II— 2.5) .  Conse¬ 
quently,  inviscid  mechanics  can  still  apply  in  Part  III  at  the  suggested  X  'v  H^  stage. 
Third,  preliminary  comparisons  with  the  computed  solutions  of  Part  I  for  small  H 
are  indicated  already.  Thus,  the  X  ^  H^  scaling  ties  in  well  with  the  computations 
(Fig.  5),  and  the  law  (II-2.4)  leads  to  P  =  H  +  o  (H)  at  eddy  closure  (where  S  -*■  0) , 
in  fair  agreement  again  with  the  calculated  results  (see  Fig.  5).  There  is  also 
agreement  concerning  the  initial  pressure  drop  (see  also  (II-2.7)  below)  and  the 
consequent  rise  initially  in  the  velocity  U(X,1)  at  the  upper  symmetry  line,  from 
(II-2.3),  as  Figs.  4-5  show.  Fourth,  since  reversed  flow  is  present,  combined  non¬ 
uniqueness  and  nonlinearity  locally  require  further  caution.  It  is  not  known  yet 
what  the  upstream  influence  of  closure  is  on  the  short  X  scale  under  discussion. 

For  instance,  a  strong  backward  jet  of  fluid  near  the  center  line  could  be  produced 
by  the  closure  downstream  (see  later  comments)  and  its  presence  when  X  is  0(A) 
would  modify  our  earlier  reasoning. 
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Nevertheless,  the  agreement  just  noted  does  give  encouragement,  albeit  guarded, 
to  a  study  of  the  description  in  (II-2.2a)  -  (II-2.5),  which  is  self-consistent  so 
far  anyway,  and  we  continue  for  now  with  the  latter.  That  description  emphasizes 
the  stage  A  =  at  which  (II-2.5)  first  becomes  significant,  leaving  the  motion 
in  sublayer  III  controlled  by  the  inviscid  system 


0  =  ii  .  p.(x) 

3Y  3X  3X  3Y  v  ' 


-  3U  - 

\p  =  t=  =  0  at  Y  =  0  [center  line], 
o  Y 


with 

1  -  -1/2 

<p  =  -  <  X  at  Y  =  S(X)  [entrainment  ] , 
Here  Y  =  HY,  and  the  expansion  implied  above, 


(II-2.6a) 


( 1 1—2 . 6b) 


3/2  -  1/2  - 

(♦,U)  =  (H  ip,  H  U) 


(H2^,  HU1) 


( I I— 2 . 6c) 


has  been  used.  Also,  the  full  symmetry  condition  is  imposed  in  (II-2.6b)  at  the 
center  line.  For  in  the  reversed  motion  there  viscous  forces  cannot  act  substan¬ 
tially  and,  in  addition,  vorticity  conservation  in  layer  I'll  coupled  with  the  initial 
zero  profile  (II-2.1d)  requires  zero  vorticity  3U/3Y  along  the  center  line  =  0+. 
Viscous  action  in  layer  III,  which  could  destroy  a  nonzero  vorticity  3U/3Y,  is 
negligible  so  far  on  this  X  ■v  scale,  being  confined  to  the  next  order  terms 
ij;^,  U]_,  which  satisfy  =  3i[^/3Y  and 


3U 

dip 

3*x 

3U  _ 

_ipi  +iis 

(II-2 ,6d) 

3X 

3X 

3Y 

3X 

3Y 

d*  3Y2 

where  H  is  the  next  order  pressure  term  in  (II-2.2a).  The  whole  local  flow 

is  then  dictated  by  the  hyperbolic  problem  (II-2.6a,  b) ,  with  (II-2.4),  the  charac¬ 
teristics  being  the  lines  ip  =  const.,  X  =  const.,  and  we  might  expect  the  condition 


i|j  ,  P  ■+  0  as  X  -*■  0. 


(II-2.6e) 


to  hold,  from  (II-2.1d). 

A  point  of  strategy  here  is  that 

P(X)  =  -  ^  U2  (X,0)  -  0 


(II-2.7) 
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from  (II-2.6a,  b,  e) ,  so  that  S(X)  >  1  from  (11-2.4)  and  wake  closure  S  -*•  0  is  not 
attainable  directly  yet,  as  remarked  before.  The  equivalent  Bernouilli  constraint 
(II-2.7)  persists  downstream  in  fact  as  long  as  the  center  line  flow  properties 
remain  predominantly  inviscid.  So  for  wake  closure  (S-0)  to  be  achieved  it  seems 
essential  for  viscosity  to  come  back  into  play  significantly  before  or  during  the 
closure  process.  Inviscid  mechanics  alone  apparently  fails  to  explain  the  eddy 
closure  here. 

The  solution  in  layer  III  still  requires  a  downstream  condition  to  be  set, 
however,  and  so  several  distinct  possibilities  arise.  These  possibilities  are 
followed  through  in  the  subsequent  sections. 
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1 1- 3.  THE  LONG  EDDY 


The  simplest  flow  solution  of  (II-2.6a,  b,  e)  is 


7  v1/2  Y 

^  -  K  X 


<  x1/2 
,  u  =  -  — 

S(X)  S(X) 


p  =  - 


1  *2 1 

(S(X))2 


(11-3.1) 


which  yields  locally  uniform  reversed  streaming.  This  motion  can  continue  for  all 
positive  X,  in  principle,  and  corresponds  to  a  downstream  condition  of  negligible 
vorticity  in  layer  III.  Allied  with  (II-2.4),  it  gives  the  implicit  equation 


S3  -  S2  =  ^  k2  X,  with  S(0)  =  1,  ( I 1—3 .2) 

for  S(X),  the  solution  being  shown  in  Fig.  7.  The  width  S(X)  of  sublayer  III  grows 
monotonically  with  distance  X  and  as  X  -*•  00 


-  P  v  S (X)  ( < 2 / 2 ) 1/3  X1/3,  0  -  -  (2k) 1/3  X1/6, 


(11-3.3) 


so  that  the  reversed  flow  becomes  accentuated.  An  eddy  closure  at  finite  X, 

O 

yielding  a  small  Xreatt  of  0(HJ),  seems  unlikely  within  such  a  structure.  Instead 
(II-3.3)  continues  downstream  until  X  becomes  0(H-^) ,  or  X  +  0(1),  at  which  stage 
the  pressure  P  ^  HP  increases  to  0(1)  and  the  outer  motion  in  layer  II  becomes 
disturbed  significantly,  from  (II-2.3). 

The  next  (and  final)  stage,  therefore,  is  relatively  far  downstream  where 
iJj,  U,  P,  X,  Y  are  all  of  0(1)  and  so  there  the  full  governing  equations  (II-2.1a-c) 
apply  across  the  entire  flow.  To  leading  order  the  starting  condition  is  now 


U-*l,  (i+YasX  +  0  for  0  <  Y  -  1 


( 1 1  —  3 . 4a) 


2  1/3  1/3 

but  with  P  'v  -  (k  /2)  X  then  and 
*  *  -  (2«,1/3  X1/6) 

U  %  -  (2k) 1/3  X1/6  j  for  0  -  Y  <  (k2/2)1/3  X1/3 


(1 1-3. 4b) 
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1/2, 


to  match  (II-3.3).  The  Chapman  singularity  is  still  present,  but  in  an  0(X  ) 

layer  centered  around  Y  =  (<^/2) ^  %!'  3  as  X  ■*  0,  and  it  requires  v'J  ^  -  <  X^^  there. 
This  is  satisfied  by  (II-3.4b).  Further  terms  in  the  small  X  expansion  (II-3.4b) 
can  be  generated  in  principle. 


The  current  stage  shows  up  a  surprising  nonuniqueness  in  the  boundary  layer 
equations,  in  fact.  For,  with  the  constraints  (II-3.4a)  holding,  one  solution  of 
(II-2.1a-c)  is  simply  the  uniform  stream  U  =  1,  \J)  =  Y  continuing  for  all  X,  Y.  So 
the  starting  form  (II-3.4b)  represents  a  second  (eigen)  solution,  branching  from 
the  incoming  uniform  stream.  Again,  the  branching  solution  appears  to  be  self- 
consistent  for  small  X  and  it  can  be  shown  to  satisfy  the  momentum  integral  equa¬ 
tion 


P(X)  + 


dY 


1 


(11-3.5) 


obtained  from  (II-2.1a),  as  well  as  the  momentum  integral  of  Part  I. 

There  appears  to  be  no  theoretical  objection  to  this  account  at  first  sight, 
nor  any  obvious  reason  why  wake  closure  and  subsequent  flow  recovery  should  not 
occur  downstream,  for  X  >  0,  yielding  self-consistency:  Fig.  7.  The  objection 
raised  earlier  by  (11-2.7)  for  instance  is  overcome  here  because  viscosity  reasserts 
its  necessary  influence  when  X  is  0(1).  Accordingly,  it  would  be  very  interesting 
to  see  the  results  of  a  numerical  study  of  the  branching  solution  stemming  from 
(II-3.4a,  b) .  On  the  other  hand,  there  is  little  evidence  yet  from  the  finite- 
spread  calculations  of  Part  I  to  support  the  prediction  of  an  0(1)  eddy  length 
as  H  -►  0.  It  could  well  be  that  the  trend  of  a  small  eddy  suggested  by  those  cal¬ 
culations  at  smallish  H  values  is  misleading  and  that  a  switch  to  the  current 
structure  with  a  longer  eddy  will  emerge  for  still  smaller  values  of  H,  but  that 
remains  to  be  seen.  Subsequently,  we  search  for  an  account  which  yields  a  smaller 
eddy,  alternative  to  that  in  (II— 3.1)  —  (II— 3.5) ,  bearing  in  mind  the  need  (satis¬ 
fied  by  the  above  solution)  for  viscous  action  during  the  closure  process  but 
noting  that  simultaneously  the  pressure  force  is  likely  to  become  negligible  then 
because  of  the  relation  (II-2.4)  with  (II-2.2a).  The  latter  property  leads  to  the 
study  in  Section  II-4  of  alternative  downstream  forms  for  sublayer  III  and  to  the 
study  in  Section  II-5  of  viscous  closure  under  zero  pressure  gradient.  Section  II-5 
has  quite  general  application  since  many  eddies  are  slender  in  practice  (e.g.,  in 
external  flow:  see  later)  and  so  cannot  produce  an  extensive  zone  of  significant 
pressure  gradient. 


*This  is  more  in  line  with  the  calculations  of  Part  1.  In  addition,  both  the 
pressure-displacement  law  (II-2.4)  at  closure  (S-*0)  and  the  momentum  integral  of 
Part  I  demand  the  same  small  pressure  rise  overall,  namely  P  =  H  +  o(H),  when  H 
is  small. 
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II-4 .  SHORTER  EDDIES 


Few  other  self-consistent  downstream  conditions  seem  immediately  possible  for 
the  sublayer  III  of  Section  II-2,  as  described  in  (II-2.6a-e)  at  least,  although 
they  do  all  tend  to  suggest  eddy  lengths  shorter  than  that  just  considered.  Two 
categories  (a)  and  (b)  arise  (see  Fig.  8),  as  described  below. 

In  the  first,  the  sublayer  III  continues  for  all  positive  X.  In  consequence, 
a  description  must  be  obtained  for  (II-2.6a,  b,  e)  when  X  ->  °°.  Since  S  ^  1  from 
(II-2.4),  (II-2.7),  there  are  only  two  plausible  trends  (labelled  (al) ,  (a2)  below), 
then:  either  S  -*■  °°  or  S  (<*>)  is  finite. 

For  subcategory  (al)  if  S  -*  »  as  X  -+  00 ,  one  solution  is  the  form  with  negligible 
vorticity,  giving  S  *  X^^,  which  was  discussed  in  the  previous  section.  Other 
candidates  arise  as  follows  (Fig.  8).  Suppose 


S  -v  s  Xn°  as  X  -*•  »,  (II-4.1) 

o 

with  the  unknown  constants  sQ,  nQ  both  positive.  Then  P  'v  -  s0  x"0  from  (II-2.4). 
Hence  (II-2.7)  implies  that  b  is  0(XnQ/,2)  sufficiently  near  the  wake  center  line, 
say  for  an  extent  0(X°1)  [n^  -  nQ]_in  the  Y  direction,  giving  ij>  of  order 
there  for  continuity,  so  that  ip  'v  xn^+n°  f(n)  with  n  =  Y/Xnl  of  0(1).  Here  f(n) 
satisfies 


n  V2  -  (n  +  2n, )  ff"  =  2n  s  ,  f(0)  =  0,  f"(0)  =  0  (II-4.2a,b,c) 

o  o  1  o  o 

+  1/2  - 

from  (II-2.6a,b).  We  may  take  it  that  f  /  -  (2sQ)  n,  for  otherwise  we  are  led 
back  to  the  earlier  solution  and  to  the  result  n0  =  n^  (=  1/3).  Again,  if  nQ  =  n^ 
it  can  be  shown  that  the  earlier  solution  is  the  only  permissible  one.  Therefore 

i  i  /  o  —  — 

we  consider  n^  <  nQ  and  f  /  -  (2s0)  n.  Then,  as  n  ->  0  (II-4.2a)  requires 
f  *  n.  So  the  derivative  of  (II-4.2a),  (nQ  -  2n^)  f*f"  =  (nQ  +  2n1)ff"’ ,  which 
gives  the  vorticity  conservation  result 


f"  -  |f|M,  (11-43) 

-M 

requires  f"  n  n  as  n  -*  0,  where  M  =  (nQ  -  2n^)/(n0  +  2n^).  In  addition,  further 
away  from  the  center  line,  as  n  -*  ®,  f  '•  -n^+  no/2ni  fr0m  (II-4.2a)  and  the  restric¬ 
tion  n-^  >  0  is  necessary.  So  ip  -Y^4n°^n^  there,  and  the  entrainment  condition 
(II-2.6b)  leaves  us  with  n0+nQ^/2n^  =  1/2,  since  S  is  0(Xn°)  from  (II— 4.1) .  Hence 
M  =  1  -  4nc,  giving  M  ^  1.  However,  f  must  be  regular  as  n  0,  to  prevent  an 
irremovable  singularity  in  higher  order  terms  in  layer  III  (see  (II  —  2.6d)  above). 


i 
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and  that  points  to  the  value  M  =  1.  Accordingly,  no  extra  solution  is  forthcoming 
here  except  perhaps  in  the  limit  as  nc  -►  0,  n^  ■  0  (n0^).  It  is  tempting  here  to 
support  the  candidate  M  =  0  (nQ  =  1/4,  n^  =  1/8),  incidentally,  associated  with  the 
simple  solution 


=  (/2  S-2  -  ft  Y/2)  Y,(l/2)ft  S2  -  S  A  S  -  2  =  <  X1/2,  P  =  1  -  S 

.1/2  ,  (,  (II-4.4) 


U  =  Jl 


t  2  K 

S  -  2  -  ft  Y,  with  S  ^  (  — 


(f) 


'v1^  V 

X  as  X  “ 


and  a  constant  vorticity  level  =  -  ft.  This  is  an  interesting  solution  which 
produces  forward  flow  along  the  center  line  and  reversed  flow  nearer  the  shear 
layer  II,  with  a  dividing  streamline  JjJ  =  0,  Y  =  2/2  S  -  2/H  returning  in  reversed 
flow  towards  the  center  line  if  ft  >  0.  Nevertheless  the  requirement  (II-2.6b)  is 
not  fully  satisfied  and,  although  the  forward  motion  along  the  center  line  can 
allow  viscosity  to  exert  a  smoothing  effect  there,  the  starting  condition  (II-2.6e) 
is  still  violated.  As  always  in  inviscid  theory,  the  possibility  of  discontinuous 
solutions  comes  to  mind,  here  and  elsewhere,  but  their  relevance  does  not  appear 
especially  likely  at  the  present  stage. 


Still  assuming  that  S  ->  «°  as  X  -*■  <*>,  and  that  the  trial  (II-4.1)  points  us  in 
the  right  direction,  we  are  left  only  with  the  question  of  the  limit  nQ  0,  then, 
with  n^  of  order  nQ^.  This  corresponds  to  a  linear  vorticity  dependence 

hz  =  X  (11-4.5) 


in  fact,  from  (II-4.3),  where  X  is  a  constant.  Thus  with  (II-4.5)  holding  another 
exact  solution  of  (II-2.6a,  b,  e)  is  obtained. 


-  -  1/2  -  -  -  1/2  1/2  - 
ip  =  -  A(X)  sinh(A  Y) ,  U  -  -  A(X)  A  cosh(A  Y) 

P  =  -  j  A  A2(X) 


(II-4.6a,b,c) 


(provided  A  >  0)  where  A(X)  is  unknown  but  A(0)  =  0.  Substituting  (II-4.6a-c) 
into  (II-2.4),  (II-2.6b),  we  have  then 


-  -  i  -2  -  -  -  1/2  -  -  -1/2 

S(X)  =  1  +  -  A  A  (X),  A(X)  sinh(A  S(X))  =  k  X  ,  (II-4.7a,b) 

two  equations  determining  S(X) ,  A(X)  for  all  X,  given  the  value  of  A  (see  Fig.  8). 

The  solution  is  regular  for  all  X  and,  as  X  ->  ”,  S(X)  grows  only  logarithmically  with 
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S(X)  ^  ~  P(X)  'v-  —  y;  +  0[ Jin  (£n  X)].  (II-4.8a,b) 

2  A1  2 

2 

A  link  with  (II-4.2a)-(II-4. 3)  for  nQ  small,  «  nQ  ,  is  then  achieved. 

The  solution  (II-4.5)-(II-4.8b)  appears  to  be  self-consistent  for  all  X.  The 
next  issue  is  whether  it  preserves  consistency  on  all  longer  scales  further  down¬ 
stream.  It  is  found  that  the  next  distinct  stage  downstream  arises  when  X  is  of 
order  H3^2(-£n  H)  .  There,  in  a  reversed  sublayer  where  Y  =  HY  and  Y  is  again 
0(1),  U  =  H1/2(-£n  H)1/2  U,  <|>  =  H3/2(-£n  H)1/2  £  from  (II-4.6a)-(II-4.8b)  and 

.  1/2  ^  ^ 

P  =  [-(-2n  H)/(4A  )  +  A  £n(-£n  H)  +  P(X)]  H  (II-4.9a) 

5/2  1/2  ^  %  1/2 

with  X  =  H  (-£n  H)  X  now,  X  >  0,  and  A^  =  1/(4A  )  is  a  constant.  The  sub¬ 

layer  therefore  achieves  an  inviscid-viscous  balance  but  is  pressure-free,  its 
governing  equations  being 


%  'v  'Xj  2^  ■v 

3U  _  dip  3U  JTU  ^ 

^  +  -2’  ^  -  ’ 
3X  3X  3Y  3Y  3Y 


and  symmetry  requires 


% 

'v  3U 

ip  =  —  =  0  at  Y  =  0, 
3Y 


(II-4.9b,c) 


(II-4.9d) 


from  (II-2.1a,b).  The  outer  condition  required  here  is  one  of  matching. 


-3/2  -3/4  1/2  -  'v  'v, 

ip  ^  -  2  A  exp[A  (Y  -  6 (X) )  ]  as  Y  >  ®  (II-4.9e) 

'V  'V 

and  needs  further  comment.  The  function  6(X)  represents  an  unknown  displacement 
effect  on  this  length  scale.  The  viscous  detached  shear  layer,  originally  layer 
I  in  Section  II-2,  is  now  centered  at  Y  =  H{ (-£nH) / 1  -A^£n(-£nH)  +  62(X)}  from 
(II-4.8a)  and  is  of  thickness  0  {H-^(-£nH)  from  Section  II-2,  where  the 
function  *s  unknown.  Between  the  detached  shear  layer  and  the  reversed 

center  line  viscous  sublayer  of  (II-4.9b-d)  the  fluid  mechanics  is  basically  inviscid 
and  pressure-free  and  so  the  ^-profile  there  must  be  purely  exponential  in  Y  to  merge 
with  layer  III  upstream.  Hence  the  form  (II-4.9e)  continues  across  the  entire  gap 
between  the  present  two  viscous  layers;  i.e.,  for  1  <<  Y  <  (-£nH)/4A^2  -Aj£n(-2nH)  + 
^(X).  At  the  upper  limit  here,  i.e.,  the  shear  layer,  under  the  assumption  that  the 
similarity  result  (II-2.5)  still  holds  in  effect,  entrainment  then  demands  that 

I 
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exp[X 


1/2 


^  % 

( <5 2 (X)  - 


^  % 

6(X))]  = 


.3/2  .3/4 


VL/2 

X 


(II-4. 10a) 


Finally,  above  the  shear  layer,  linearized  inviscid  channel  flow  theory  applies 
again,  giving  the  relation 


P(X)  =  -  62(X).  (II-4 . 10b) 


The  rationale  here  is  that  the  viscous  sublayer  problem  (II-4.9b-e)  is  the  central 
one,  determining  <5^  hence^(II-4 . 10a)  determines  &2>  then  (II-4.10b)  determines  the 
pressure  response  P,  and  P  then  drives  higher  order  corrections  to  the  flowfield. 

a, 

The  rationale  works  for  small  positive  X,  where  a  match  with  the  earlier 
properties  for  X  finite  is  achieved.  Thus  the  reversed  sublayer  (II-4.9b-e)  be¬ 
comes  predominantly  inviscid  then,  and  is  regular,  with 


*  ^  F  (Y)  +  X  F  (Y)  +  0(X2)  (11-4. 10c) 

o  1 

where  Fq  =  -2  X  sinh  (X^2  Y)  from  (II-4.6a)-(II-4.8a) .  Here  (11-4. 9b-d) 
requires  that 

1/2  -  -1  1/2  - 

F^(Y)  «  cosh(X  Y)(tan  [exp(X  Y)  ]  -  tt/4)  (II-4.10d) 

and  so  we  have 


i^i  'V  'V 

|  6  |  «  X  as  X  -*■  0  +. 

Further  terms  can  be  generated  for  small  X. 


(II-4 . lOe) 


The  description  apparently  cannot  continue  for  all  positive  X  however.  A  sim¬ 
ilarity  solution  does  suggest  itself  for  X  -*•  °°,  in  the  form 


£  -v  X  F(Y)  +  0(1)  (II. 4-11) 

% 

which  yields  for  F  the  nonlinear  differential  equation 

F"’  +  F  F"  -  F’2  =  0  with  F(0)  =  F"  (0)  =  0  (II-4.12a-c) 

from  (II-4.9b-d).  But  an  integral  of  (II-4.12a-c)  is 
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F"  =  exp 


Y  Y  i 

[ -  J  F(yx)  dy1]  J  F'2(y2)  exp[ J  F(y1>  dy1]  dy2  (II-4.13a) 


so  that  F"  >  0  for  Y  >  0.  Hence  the  required  asymptote. 


F  'V/  -  2  X  exp[X^2  (Y  -  const.)]  as  Y  -*•  “  (TI~4.13b) 

a. 

from  (II-4.9e),  is  unobtainable  since  (II-4.13b)  requires  F"  <  0.  Instead 
(II-4.12a)-(II-4.13a)  can  be  shown  to  give 

-v  1/2  - 

F  'v  a  exp(X  Y) ,  Y  -*■<*>,  a^  >  0  •  (II-4.13c) 

In  a  sense  it  is  a  pitv  that  (II-4.13b)  cannot  be  satisfied  with  (II-'i .  1 2n-c) .  For 

'V  -1/9  9, 

(I 1-4 . 13b) ,  with  (11-4,11),  implies  that  6  %  -X  LU  In  X  as  X  +  »;  so  (II-4.10a,b) 

r\j  y  1/2  % 

give  -P  =  ^2  ^  (1/2  A  )£n  X  ->  -°°  and  the  detached  shear  layer  approaches  the 
center  line  at  last;  the  shear  layer  and  the  reversed  center  line  sublayer  then 
merge  (anticipating  closure)  downstream  when  X  is  0(H2);  the  boundary  layer  equations 
apply  there  without  pressure  gradient;  and  self-consistency  could  be  achieved.  The 
nonexistence  of  a  solution  to  (II-4.12a-c)  with(II-4 . 13b)  rules  out  such  an  account, 
of  course,  but  it  does  raise  again  the  question  of  wake  closure  occurring  without 
significant  pressure  forces:  see  Section  11-5  below. 

For  the  reversed  sublayer  of  (II-4.9b-e)  the  solution  most  likely  terminates  in 
a  singularity  at  a  finite  positive  value  Xn  of  X.  The  singularity  as  X  +  X  -  is  given 
by  (II-4.11)  again,  effectively,  with  X  replaced  by  (X-XQ)  and  self-consistency  is 
maintained  now  with  (II-4.13c)  yielding 

•v  'v  -1/2  <v  'v  'v  'v 

6 (X)  %  -X  £n  (XQ  -  X)  as  X  ■+  XQ-.  (II-4.14) 

Hence,  6  and  £2  +  00  then,  the  detached  shear  layer  departs  further  from  the 

center  line  and  the  pressure  response  ^  -*■  -°°,  from  (II-4. 10a, b) .  The  further  pres¬ 
sure  fall  due  to  this  departure  from  the  center  line  eventually  affects  the  flow 
in  the  reversed  sublayer  first,  within  a  more  local  scale  around  X  =  XQ,  and  so 
then  the  relations  (II-4.9e),  (II-4.10a)  become  altered.  No  clear  means  of  closure 
of  the  eddy  in  these  circumstances  has  emerged  yet,  but  the  pressure  response  could 
cause  the  detached  shear  layer  to  reverse  its  trend  abruptly  and  turn  back  towards 
the  center  line  within  the  localized  scale.  For  no  solution  of  the  sublayer 

'V  'V 

(II-4.9b-e)  exists  just  downstream  of  the  singular  point  X  =  XD,  and  so  a  change  of 
structure  is  necessary  locally.  If  that  causes  wake  closure  the  prediction  Xreatt  = 

I  0[H5/2(-HnH)1/2]  follows. 

I 

1  So  far,  therefore,  only  two  accounts  have  avoided  obvious  failure,  the  one  of 

Section  1 1  —  3  and  the  one  just  described  in  ( 1 1-4 . 5)- ( 1 1— 4 . 14) . 
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A  third  candidate  emerges,  subcategory  (a2),  when  S(°°)  is  assumed  to  be  finite. 
This  yields  the  property 


ip  ->  \p  (Y),  U  ->  U  (Y),  P  ->  1  -  S(°°)  as  X  -*  »  (II-4.15a) 

oo  oo 


in  layer  III  of  Section  II-2.  Here  the  profiles  ^(Y) ,  U^Y)  =  (Y)  are  arbitrary 

except  that  1^(0)  =  ^"(0)  =  0  for  symmetry,  we  might  expect  U^CO)  <  0  and,  because 
of  the  increasing  entrainment  (II-2.6b),  ^(Y)  +  -®  as  Y  +  S(°°)  So  the 

reversed  velocity  in  III  becomes  accentuated  just  below  the  shear  layer  I. 

The  nature  of  the  reversed  acceleration  here  controls  what  happens  further 
downstream  on  longer  length  scales.  Suppose  that 


(Y)  -b  (S<«)-  Y  _Yas  Y  -*•  S(“)  - 

oo  X 

where  b^,  '<  are  positive  constants.  Then  as  X  -*•  °° 


(II-4 . 15b) 


1  -  S  -  P  *1  -  S(°°)  +  b„  X  1/2',  ip  ^  ip  (Y)  +  x'1/2>  J  (Y)  (II-4 . 15c) 

l  00  A 


and,  from  (II-2 . 6a,b ,e) , 


VY)  =  -b2 


x 

5-«>  /  o% 

OO 


( I 1—4 . 15d) 


with  the  constant  b2  remaining  unknown.  Hence  ^  induces  an  increasing  displace¬ 
ment  effect  as  Y  -*•  §(<*>)-.  As  a  result  (Fig.  8),  in  a  thinner  layer  just  below 
Y  =  S(®)-,  the  solution  is 


-1/2  -  -7 

<1/  ^  -b  X  [b  I  -  n]  (II-4.15e) 

1  2  A 

F(°°) 

for  n  =  (Y  -  5 (-))  X  (  now  of  the  order  unity.  Here  the  constant  IA  (Y)dY 

is  finite  and  positive,  while  (II-4.15e)  merges  with  (II-4.15c,d)  as  n  ■*  and 
enforces  the  entrainment  condition 


(1 


+  iAr  =  bY  < 


( 1 1—4 . 1 5  f ) 
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from  (II-2.6b),  since  the  shear  layer  is  at  n  =  -^>2 •  Equation  (II-4.15f)  fixes  the 
value  of  b2  in  terms  of  the  final  profile  ^(Y),  since  the  latter  determines  b^  and 
1^.  Further  contributions  can  be  generated,  in  inverse  powers  of  X,  and  forms  other 
than  (II-4.15b)  are  also  possible,  but  the  main  point  is  that  again  a  self-consistent 
description  of  layer  III  emerges. 

Can  this  third  account  also  remain  consistent  further  downstream,  producing 
wake  closure?  The  next  distinct  stage  arises  when,  in  layer  III,  the  maximum 
(reversed)  velocity  rises  to  0(1)  because  then  the  forward  shear  layer  I  must  be 
affected  significantly  by  the  motion  underneath  it.  This  maximum  occurs  where 
(II-4.15e)  holds,  giving  U  of  order  /2y ,  i.e.,  U  of  order  x(y+1)/2y  f rom 

(II-2.6c),  when  X  >>  1.  Consequently,  since  X  =  X,  the  stage  X  =  0[H^  »+3)/(y+  )]^ 

o  9 

greater  than  H  but  less  than  H  ,  is  the  next  crucial  one  downstream.  Note  that  the 
thin  layer  where  (II-4.15e)  holds  is  of  thickness  0[H  x“1/2Y]  ,  i.e.  ,  0[H^+^2Yx-1/2y j . 
hence  at  the  stage  just  mentioned  the  thickness  decreases  to  0 [H^Y+3)  ^  (2y+2)  ^ , 
exactly  the  viscous  thickness  of  0(X^^)  associated  with  the  Chapman  solution. 

The  latter  solution  can  apply  no  longer  then  and  the  whole  viscous  shear  layer  must 
change  in  character.  This  prospect  is  studied  further  in  Section  II-8  below. 

The  second  category,  category  b,  has  the  sublayer  III  terminating  after  a  finite 
distance,  at  X  =  X  -  say,  suggesting  Xreatt  =  O(H^).  Terminal  forms  can  be  written 
down  which  provide  a  self-consistent  description  of  layer  III  as  X-+XQ-.  The 
mechanism  for  completing  eddy  closure  in  such  cases  remains  elusive,  however,  mainly 
because  of  the  objection  (to  inviscid  theory)  arising  from  (II-2.7).  Viscous  forces 
must_exert  a  strong  influence  before  closure  can  take  place;  yet  a  termination  at 
X  *  XQ  provides  no  clear  sign  of  increased  viscous  action  there.  In  contrast,  the 
arguments  supporting  the  earlier  possibilities  of  Section  II-3  and  (al) ,  (a2)  above 
all  bring  viscosity  back  into  play  before  closure.  A  number  of  alternatives  arise 
here  including:  the  different  forms  of  viscous  shear  layers  (Section  II-8) ;  higher 
order  terms  in  III,  e.g.,  as  in  (II-2.6d),  may  reinstate  viscous  effects  near  the 
center  line  as  X-*XQ  if  P,  U-K)  there;  layer  III  may  have  a  discontinuous  inviscid 
solution;  multiple  eddies  may  occur,  e.g.,  if  X  in  (II-4.5)  is  negative ;  a  backward 
jet  may  be  present  along  the  center  line.  None  of  these  has  been  found  to  produce 
a  complete  self-consistent  argument  for  a  finite  X  termination  yet,  although  Section 
II-8  below,  concerning  nonentraining  shear  layers,  does  revitalize  the  possibility. 
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I I— 5 .  ON  PRESSURE-FREE  VISCOUS  WAKES 


Viscosity  has  to  come  back  into  the  reckoning,  for  wake  closure  to  take  place, 
but  if  in  addition  the  length  scales  governing  the  closure  are  relatively  small 
(unlike  Section  II-3)  then  the  induced  pressure  is  small,  in  view  of  (II-2.4). 
Negligible  downstream  pressure  is  an  attribute  of  the  new  inviscid  solutions  of 
Section  II-4,  for  instance.  Since  the  shear  layer  velocities  are  of  0(1)  by  con¬ 
trast,  the  properties  of  pressure-free  viscous  closure  must  be  addressed. 

Our  reasoning  here  hinges  on  the  flow  features  near  the  centerline,  given  the 
pressure-free  equations  (for  >  0) 


2 

9y  U3U  _  9<li_  _9U _ +  9  U 

9y  ’  9x,  9x  9y  .  2 

1  t  11  ay^ 


(11-5. la) 


holding  on  some  short  length  scales  X  =  Aj  x^  ,  y  =  Aj  y^,  with  U  =  0(1).  Here 
1  is  unknown  as  yet,  but  A-.  >*  A  and  the  constraints  are 


t P  =  —  -  0,  y  =0  (symmetry),  and  U  *  1 

dyl  1 

as  y^  "  (freestream) . 


(11-5.  lb, c) 


The  constraints  stem  from  (II-2.6b)  and  (11-2.3),  respectively,  for  the  wide-spread 
cascade,  but  the  application  intended  is  broader  since  many  other  wake  flows  with 
closure  seem  too  slender  to  provoke  a  substantial  pressure  force. 

If  the  starting  profile  at  x-^  =  0+  has  only  forward  velocities  then  (II-5.1a-c) 
can  be  integrated  forward  in  x^  ,  giving  a  standard  wake  problem.  If  i losure  has 
still  to  occur,  however,  the  starting  profile  must  contain  some  reversed  velocities. 
We  consider  whether  such  reversed  flow  can  be  present  near  the  center  line  (to  set 
up  a  downstream  condition  for  a  layer  such  as  III  above)  and  be  initially  small  (to 
avoid  sending  a  strong  reversed  jet  back  into  the  earlier,  shorter  length  scale, 
motions).  That  indicates  a  local  similarity  solution  holding  as  xj  •0+.  of  the  form 
tfi  =  /(K+1)G(0  where  '  =  Here  K  is  an  unknown  constant  but  K’  1 

so  that  U  =  :'i (K-l)/ (K+1)q» (r)  is  small.  If  K  1  this  viscous  region  is  contracting 
as  xjy>-0  ,  whereas  it  is  expanding  if  K  -  -l.  The  function  C  satisfies 


j(K  +  1)  G'"  +  K  G  G"  -  (K  -  1)  C 


G (0 )  =  G"(0)  =  n 


(I  1-5. 2a) 


(I  1-5. 2b) 
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from  (II-5.1a,b).  An  integral  of  (II-5.2a,b) 

G"  "  fjl  eXP  Si  G'2«2>  “P  |Sl 

(11-5. 3a) 

then  shows  that  the  vorticitv  is  positive, 

G"  >  0  for  (  >  0,  (II-5 . 3b) 

since  |k|  >  1.  Therefore  the  velocity  G' (£)  increases,  with  £,  from  the  unknown 
center  line  value  G'(0).  As  £  continues  to  increase  only  two  ultimate  trends  can 
be  observed:  see  Fig.  9.  For  K  <  -1  a  singularity  is  encountered  at  a  finite 
position  £  =  £Q  >  0,  with 

—  I K I  —1 

G  O,  +  C  (£  -  O  1  '  +  0(5  -  £)  as  £  -  £  ■  (II-5.4a) 

1  o  o  o’ 

whereas  for  K  >  1  the  solution  continues  for  all  £,  giving  (apart  from  an  origin 
shift) 


G  +  C2 


K 

€  + 


0(C_1)  as  £ 


(II-5.4b) 


Here  and  C2  are  unknown  constants  but  both  are  positive  because  of  (II-5.3b),  so 
that  the  ultimate  trends  produce  forward  flow,  G'  -»■  +  <*>.  It  is  interesting  that  the 
forms  (II-5.4a,b)  are  inevitable  for  any  starting  value  G’(0),  whether  positive  or 
negative  (the  case  G'(0)  =  0  gives  the  trivial  result  G  =  0  and  is  of  no  concern). 
This  implies  the  surprising  result  that ,  for  instance ,  the  near-wake  solution  of 
Goldstein  (Ref.  11)  corresponding  to  K  =  2  above  is  not  unique  (see  Ref.  17):  if 
reversed  center  line  flow  G' (0)  <  0  is  allowed  then  a  second  solution  exists  which, 
like  Goldstein's,  still  gives  G"(<*>)  finite  and  positive. 

Suppose  first  that  K  <  -1,  the  expanding  case,  implying  the  property  (Ii-^.4a). 
A  slimmer  region  of  thickness  0(1)  in  y^  is  induced  (see  Fig.  9)  near  f,  =  £0  to 
smooth  out  the  growth  of  (II-5.4a).  The  solution  there  has  faster  flow  with  U, 
of  0(1) ,  say 
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(o)  —  + 

ip  ip  (y)  as  0 


(II-5 . 5a) 


where  y-^  =  x-^  ^  +  y  is  large.  The  governing  equations  (II-5.1a)  are 

satisfied,  and  the  profile  (y)  remains  arbitrary  apart  from  the  requirement 


(y)  ^  C2  I  y  j  ^  as  y  -» 


( II— 5 . 5b) 


of  matching  with  (II-5.4a),  at  one  extreme,  and,  at  the  other,  a  matching  condition 
required  as  y  -*■  00 .  The  latter  depends  on  how  multi-structured  the  local  wake  is. 
The  simplest  condition  is 


^^°\oo)  =  1 


(II-5.5c) 


which  satisfies  (II-5.1c)  directly  and  gives,  for  the  wake  displacement  6(x^) 
defined  by  ip  'v  y^  -  6(x^)  in  (II-5.1a-c),  the  growth 


6(x,)  'v  C  x 

o  1 


-1/ ( I K I -1) 


=)  as  x 


(11-5. 5d) 


This  growth  upstream  appears  physically  sensible  with  regard  to  both  the  wake 
closure  anticipated  further  downstream  and  to  the  detached  flow  properties  (Sections 
11-2,4)  holding  on  shorter  scales.  So  the  question  arises  of  whether  those  shorter 
scale  flows  can  be  joined  to  the  present  viscous  structure  for  small  x^,  or  not, 
given  that  in  particular  the  0(1)  thick  region  of  fast  flow  in  which  (II-5.5a) 
holds  must  exhibit  the  algebraic  decay  (II-5.5b)  with  | K |  >  1.  Again,  realistic 
conditions  other  than  (II-5.5c)  are  possible  for  y  -*•  <*>.  For  instance,  a  decay 
like  that  in  (II-5.5b)  can  be  achieved  for  y  ->•  °°  with 


’/°\y)  ^  C3  +  C4  y  v 


( 1 1— 5 . 5e) 


where  v  (>  0),  and  (a  mass  flux)  are  constants.  If  (II-5.5e)  holds  then  a 
further  thicker  viscous  region  arises  above  the  fast  flow  region  with  thickness 
y  =  0(x^  ^  ^)  and  ip  -  C3  =  0(x^~^''  ^  )  is  small  provided  v  >  1.  In  effect, 

the  governing  equation  then  is  (II-5.2a)  again,  with  K  replaced  by  -v  and  ( 1 1— 5 . 2b ) 
is  replaced  by 


G  -  c4  a  -  e0) 


(11 . 5— 5  f ) 
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from  (II-5.5e).  As  £  increases  further  the  solution  tends  to  one  of  two  forms, 
both  yielding  slower  motion  (see  Ref.  17,  however), 

G  'v  6  £  ^  or  const.,  as  £  -*•  °°,  (II-5.5g) 

consistent  with  (II-5.2a);  or  a  singularity  like  (II-5.4a)  is  repeated.  Conversely, 
if  v  <  1  then  no  such  outer  region  arises.  This  enables  the  present  form  to  be 
joined  directly  to  a  further  incoming  shear  layer  more  removed  from  the  center  line 
if  the  mass  flux  ^  approaches  the  value  in  the  lower  reaches  of  the  shear  layer. 
There  is  a  strong  connection  between  these  properties,  involving  algebraic  decay 
at  the  onset  of  wake  closure,  and  those  of  Section  II-4(a2),  where  algebraic  decay 
was  found  to  be  permissible  nearer  the  start  of  the  motion,  and  Section  II-8  below 
pursues  that  theme. 

If  K  >  1,  on  the  other  hand,  the  relevance  to  any  shorter  scale  detached  flow 
appears  less,  in  view  of  the  feature  (II-5.4b)  where  C2  >  0  gives  forward  motion. 

Thus  the  majority  of  the  wake  is  initially  forward  on  the  present  scale  and  major 
adjustments  of  the  earlier  detached  mption  would  need  to  have  taken  place  already, 
on  a  shorter  scale.  In  particular,  a  join  from  the  center  line  y^  =  0  to  an  incoming 
detached  shear  layer  is  unlikely  with  K  >  1  unless  the  lower  part  of  the  shear 
layer  has  acquired  substantial  forward  motion  beforehand. 

Other  issues  arising  here  are  as  follows.  First,  the  limit  as  K  +  ■»  above 
recovers  the  finding  of  Section  II-4(al),  in  (II-4.11)-(II-4.13c) ,  where  the  equa¬ 
tion  is  related  to  that  of  the  asymptotic  suction  profile.  Second,  simple  power 
solutions  for  G(£)  in  (II-5.2a)  exist  only  for  K  =  1,2.  For  K  =  1, 

G(£)  =  *  const.  £  satisfies  all  of  (II-5.2a,b)  but  this  corresponds  to  U  =  *  const., 
producing  a  centerline  flow  which  is  not  coming  to  rest,  in  contradiction  to  the 
shorter  scale  features.  For  K  =  2,  G(£)  =  -  d^£^  satisfies  (II-5.2a)  but  not 
(II-5.2b);  instead,  it  satisfies  a  no  slip  condition,  a  matter  taken  up  in  Section 
II-6.  Third,  although  (in  consequence)  a  study  of  the  limit  K-+1+  suggests  itself, 
further  analysis  shows  that  still,  with  virtually  no  pressure  gradient  present,  a 
pronounced  interval  of  reversed  flow  remains  impossible,  in  line  with  (II-5.4b). 

If  a  significant  pressure  gradient  is  acting  then  more  possibilities  open  up  but 
really  that  forces  us  straight  back  to  the  self-consistent  solution  of  Section  II-3 
and  its  scalings,  in  contrast  with  the  smaller  scale  phenomena  being  sought  here. 
Fourth,  quite  different  self-consistent  accounts  for  the  start  of  the  viscour  wake 
(II-5.1a-c)  are  obtainable,  with  the  centerline  motion  reversed  as  ^|^uired,  if  U,  ip 
are  taken  to  have  finite  nonzero  profiles  for  y^  of  0(1),  say  p  -*■  <J>  ly^),  as 

x^  ■*  0  .  This  is  similar  to  (II-5.5a).  The  structure  on  a  shorter  length  scale, 
however,  then  has  to  contend  with  a  strong  backward  jet  at  the  centerline,  and  that 
is  a  major  task:  see  Sections  II-7,  II-8  below. 
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Next,  an  integral  of  (II-5.1a-c)  requires  the  momentum  deficit  integral  to 
be  conserved, 


U(1  -  U)  dy: 


o 


I 

o 


(II-5.6a) 


where  I  is  a  constant,  deducible  from  the  starting  form  of  the  motion.  Far 
downstream  if  wake  recovery  occurs,  then  the  left  hand  side  of  (II-5.6a)  tends  to 
6  (°°) ,  so  that  6  (°°)  =  I0.  But  the  linearized  relation  (II-2.4)  still  applies  in 
effect  here,  i.e.,  P  =  H  -  6(x^)  for  all  x-^  since  the  present  length  scales 

are  short.  Hence 


P  -*•  H  -  I  as  x, 

o  1 


(11-5. 6b) 


The  corresponding  integral  of  the  entire  problem  (II-2.1a-d)  requires 


2  3 

P  -►  H  +  H  +  H  +  .. 


(II-5.Sc) 


>  2 

far  downstream,  however.  Therefore  two  cases  emerge.  Either  -  0(H  )  (yielding 
Xreatt  -  in  which  case  the  starting  solution  for  (II-5.1a-c)  must  satisfy 

IQ  =  0,  which  is  a  severe  constraint.  Or  A^  <<  ,  ^reatt  <<r  ,  in  which  case  I0 

can  be  nonzero  if  A-^  =  H^  but  otherwise  must  be  zero  again.  Notice  that  IQ  cannot 
be  positive,  from  (II-5.6b,c).  This  singles  out  the  scaling  A^  =  n  ,  ^reatt  = 
if  the  integrated  momentum  deficit  IQ  is  initially  nonzero  and  negative.  Moreover, 
if  IQ  <  0  then  initially  there  must  be  substantial  regions  of  strong  reversed  flow 
(U  <  0)  or  of  velocity  overshoot  (U  >  1),  because  of  (II-5.6a).  In  any  region 
where  |u|  is  small  the  integral  in  (II-5.6a)  contributes  only  the  jump  in  tp  across 
the  region  and  this  is  negligible  from  the  starting  forms  of  (II-5.2a-c).  Some 
care  is  necessary  over  these  conclusions  when  6  (o+)  is  not  finite,  as  in  (II-5.5d), 
of  course,  since  ( 1 1—5 . 6b )  becomes  questionable  and  the  connection  back  to  the 
earlier  shorter  scale  properties  then  is  not  yet  clear;  nevertheless,  the  xesults 
(II-5.6b,c)  must  always  be  reconciled. 

The  sixth  and  final  issue  here  concerns  the  integral  form 
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obtained  from  (II-5.1a-c),  where  V  =  -  3^/Sx^:  c.f.,  (II-5.3a).  If  the  centerline 

flow  is  coming  to  rest  at  -*•  0+,  U(xpO) -*•  0,  as  supposed  prior  to  (II-5.2a),  then 
(II-5.7a)  yields  3U/3yj  >  0  for  small  y^  >  0.  So  the  velocity  U  has  a  local  minimum 
at  the  centerline. 
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II-6.  ON  MASSIVE  EDDY  CLOSURE  AND  REATTACHMENT 
ONTO  A  SOLID  SURFACE 


This  section,  a  by-product  of  Section  II-5,  hinges  on  the  feature  that  the 
exact  solution  of  (II-5.2a)  with  reversed  flow  for  K.  =  1  is  G  =  -d^  giving 


■-2diV 


"‘"Vi 


(II-6.1a) 


the  constant  d^  is  assumed  positive  here.  The  local  solution  (II-6.1a)  for  small 
Xi  satisfies  the  following  conditions:  (1)  no  slip  holds  at  y-^  =  0;  (2)  the 
Chapman  entrainment  ip  =  -k  x^  is  supplied  at  the  position 

y^  =  (K/d^^)1^2  x  ^  + - -  ( 1 1—6  .  lb) 

1/2 

around  which  the  thinner  0(x^  )  Chapman  form  applies;  (3)  the  reversed  velocities 

U  induced  between  the  centerline  y^  =  0  and  the  shear  layer  at  (II-6.1b)  are  small, 
of  order  xj^^  as  x-^-K),  in  keeping  with  an  assumption  of  slower  flow  beneath  the 
shear  layer  on  shorter  length  scales;  (4)  the  pressure  response  required  remains 
negligible  (see  comments  on  previous  section) . 

Hence  the  initial  form  (II-6.1a,  b)  of  the  boundary  layer  (II-5.1a,  c)  can 
describe  the  start  of  viscous,  pressure-free,  reattachment  onto  a  solid  surface, 
in  principle,  with  (II-5.1b)  replaced  by  ^  =  U  =  0,  y^  =  0.  See  Fig.  10,  also 
Ref.  17.  Further  terms  in  the  small-x^  expansion  stemming  from  (II-6.1a,  b)  can 
be  generated  at  will,  including  exponentials  (see  Section  II-7  and  Appendix  D 
below),  integral  properties  like  those  of  Section  II-5  can  be  derived,  and  it 
would  be  interesting  to  follow  the  subsequent  downstream  development  of  (II-5.1a,  c) 
with  no  slip  at  y^  =  0,  which  necessarily  poses  a  challenging  numerical  task. 
Meanwhile,  we  consider  below  the  repercussions  of  (II-6.1a,b) ,  and  its  features 
(l)-(4)  above,  for  the  shorter  scale  motion.  The  applications  in  mind  here  include 
the  gross  separation  and  subsequent  reattachment  in  flow  past  ramps,  past  bluff 
bodies  with  trailing  splitter  plates,  and  in  wind-tunnels,  all  of  which  are  of 
wide  concern. 

The  first  distinct  stage  encountered  as  the  length  scale  is  shortened  (say 
x^->0(A2),  small)  appears  where  the  induced  pressure,  while  still  small,  first 
affects  the  centerline  motion,  since  the  latter  is  relatively  slow  and  decelerating, 
from  (II-6.1a,  b) .  There  the  boundary  layer  equations  still  hold,  with  induced 
pressure  included,  in  a  thin  sublayer  where  y^  is  0(^2  and  U,  ip,  P  become 
0(A21/3,  a22/3,  a22/3)  in  turn.  So,  in  effect,  Eqs.  (II-2.1a)  apply  again  but  now. 


UTRC83-13 


U  =  \p  =  0  at  Y  =  0  (no  slip) 


( 1 1-6 . 2a) 


^  'v  -d^  Y  as  Y  +  “  (outer  match) 

2  2/3 

ip  ^  - d 5  X  as  X  +  +  <»  (downstream  match) 


where  £  -  Y/X  ,  and  the  induced  pressure  is 


(11-6. 2b) 


(11-6. 2c) 


p(x)  =  -  (x/d^172  x174, 


( 1 1-6 . 2d) 


Here  (see  Fig.  10)  the  pressure  is  provoked  by  a  combination  of  the  interaction 
law  P  =  -S  (due  to  the  current  small-scale  effect  on  the  outer  flow,  as  in 
(II-2.4)),  and  the  entrainment  law  ip  -  -irX^7^  as  Y  +  S  -,  from  (II-2.5),  giving 
d^  S2  =  k  X172  because  of  (II-6.2b).  This  combination  is  analogous  to  that 
in  Section  II-4  (al)  and  is  taken  only  as  an  example,  relating  to  the  cascade  with 
a  splitter  plate  present  along  y  =  0.  The  other  applications  mentioned  above 
yield  pressure  forms  different  from  (II-6.2d),  although  equally  simple  in  certain 


The  boundary  layer  problem  (II-6.2a-d)  with  (II-2.1a)  looks  unusual  at  first 
sight,  in  having  its  starting  condition  at  downstream  infinity,  in  the  uniform 
shear  flow  (II-6.2c),  but  that  is  because  the  motion  is  reversed.  The  flow 
solution  can  develop  satisfactorily  there  in  the  form 


ip  ^  -dj'V73  +  Xl74gl(C)  + - -  X 


( I I — 6 . 3a) 


Here  g^  satisfies  g^'''  -  d^  (2C2g1',/3  -  >,  gj/2  +  g-^/2)  =  -(x/lbd^) 
gl(0)  =  g{(°)  =  gj'C00)  =  °»  from  (II-2 .  la) ,  (II-6.2a,  b,  d) .  Hence 


1/2  with 


gj'  =  (27x) 172  J  exp(-rd1^ 3) (2+9r) 


-3/4  -7/12  ,  ,fn,  5/6 


dr/  [  12 (-1/ 3)  .'d, 


(11-6. 3b) 


and  g^'  (0)  ;•  0.  The  incident  reversed  flow  downstream  is  therefore  decelerated 
further  by  the  pressure  field  as  X  decreases.  Further  terms,  including  an  origin 
shift  in  X,  can  be  generated  in  the  expression  (II-6.3a). 


Since  the  local  problem  (II-6.2a-d),  (II-2.1a)  has  a  prescribed  pressure 
gradient  which  is'adverse  in  the  minus  X  direction,  however,  the  flow  solution 
cannot  continue  for  all  X  and  come  to  rest  at  X  =  0  +  as  the  body  scale  properties 
demand.  Instead,  either  the  classical  Goldstein  (Ref.  12)  separation  singularity 
is  produced  at  a  finite  positive  value  of  X,  say  X-Oi^  +  ,  or,  due  to  an  origin 
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shift  if  allowed,  the  flow  reaches  X  =  0  +  with  nonzero  velocity.  We  expect  the 
former  option  in  general.  Fortunately,  Goldstein’s  singularity  is  not  necessarily 
a  disastrous  occurrence  here,  if  around  X  =  local  interaction  brings  in  the 
pressure-displacement  law  (2.4)  which  yields  regular  separation,  from  Ref.  5. 
Therefore  secondary  separation  probably  can  be  accommodated  either  directly  or 
interactively,  depending  on  the  particular  local  interaction  near  X  =  X-^, 
and  self-consistency  seems  not  unlikely.  Moreover,  in  some  applications  the 
pressure  form  replacing  (II-6.2d)  turns  out  to  have  favorable  gradient,  thus 
preventing  secondary  separation,  i,i  aspect  worth  following  through. 

The  overall  trend  of  the  flow  structure  produced  above  as  a  result  of 
(II-6.1a,  b)  is  encouraging,  then,  as  a  strong  back  flow  is  prevented  by  means 
of  viscous  action  but  little  pressure  force.  There  seems  good  reason  to  pursue 
the  matter  further,  relating  to  the  specific  applications  mentioned  before,  but 
that  must  be  regarded  as  future  work.  We  return  now  to  the  original  wake  closure 
problem. 
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I 1-7.  CONNECTION  FROM  WIDESPREAD  CASCADE  TO  EXTERNAL  FLOW; 
AND  HOW  FAR  CAN  THE  KIRCHHOFF  PARABOLA  LAST? 


Before  considering  the  relation  between  the  widespread  cascade  and  external 
flow,  we  discuss  certain  features  of  the  external  flow  problem  itself. 

In  external  flow  past  a  bluff  body,  Kirchhoff’s  inviscid  separated  flow 
solution  produces  a  parabolic  growth  (Fig.  11) 


S(x)  ^  bx 


1/2 


as  x  -►  °° 


(II— 7 . 1) 


for  the  eddy  height  y  =  S(x)  downstream,  as  mentioned  in  Part  I.  The  viscous  shear 

—1  /  2  ~ 

layer  near  the  eddy  boundary,  specifically  where  y  =  S(x)  +  Re  y  with  y  of  0(1), 
also  grows  parabolically  as  x  +  ®  if  the  Chapman  form  (II-2.2d)  holds  then,  with 

n  =  y/x1/2,  u-fc'(n),  *  *  Re"1/2  x1/2  fc(n),  fc'(«)  =  1,  fc’(-“)  =  0,  ft(— )- 
Also,  Ref.  1  shows  that  the  mild  uniform  reversed  streaming 


u 


Re 


-1/2 


b 


,  „  -1/2  Ky 
^  ^  -Re  -g— 

0  £  y  <  bx1/2. 


as  x 


( I 1—7 .2) 


then  satisfies  the  inviscid,  but  slower  flow,  properties  of  the  eddy  as  well  as  the 
entrainment  (^  -*■  -k  Re-l/2  x1/2,  y  -*•  bxl/2  -  )  and  symmetry  (y  ■+  0)  conditions, 
downstream. 

Viewed  separately,  each  of  the  three  parts  of  the  downstream  flow  here,  (1)  the 
uniform  stream  u  =  1  above  the  shear  layer,  (2)  the  Chapman  form  within,  and  (3) 
the  mild  reversed  stream  (II-7.2)  underneath,  constitutes  an  acceptable  downstream  limit 
for  the  Navier-Stokes  equations,  and  so  produces  no  call  for  wake  closure  to  occur  at 
all.  Yet,  viewed  together,  (1)— (3)  do  not  satisfy  the  Navier-Stokes  equations  in  a 
far-wake  sense,  for  there  (II-5.2a,b)  apply  with  K  =  1,  so  that  (II-5.3a)  then  requires 
zero  vorticity  G"  and  only  the  limit  u  =  1  remains  acceptable.  Therefore  wake  closure 
has  to  occur,  and  the  (open  wake)  description  involving  (1)— (3) ,  stemming  from 
(II-7.1),  must  break  down  on  some  long  length  scale  downstream.  It  cannot  persist 
indefinitely. 

This  brings  us  to  the  question  posed  above;  what  is  that  length  scale?  Also, 
what  forces  the  breakdown  to  occur?  The  former  question  is  important  when  the 
relation  between  the  widespread  cascade  and  external  flows  is  considered,  while  the 
latter  question  leads  to  some  perhaps  unexpected  aspects  of  separated  motions. 

The  reasoning  involved  must  be  tentative  again,  because  recirculatory  motion  is 
present,  but  consistency  overall  is  the  eventual  aim. 
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The  matter  centers  on  certain  exponentially  small  terms.  Take  the  body  scale 
flow  first  (see  Fig.  11).  In  the  eddy  (E),  an  acceptable  solution  is  (iji,  u,  p)  ■ 
(Re“^^5,  Re-^/^  Re“l  p  +  const.)  to  leading  order,  with  x,y  of  0(1),  so  that 

u,  p  are  governed  by  the  Euler  equations.  Vorticity  conservation  gives  V2  vp  =  f(^) 
and  if  (II-7.2)  applies  then  f  =  0.  Also, 


=  ~<(x) 


at  y  =  S(x)-,  ^  =  0  at  y  =  0 

°y 


(11-7. 3a) 


from  entrainment  and  symmetry  respectively  y(see  also  Section  II-2).  Here 
K(x)  s-4i(x,-co)  is  the  reduced  mass  flux  required  by  the  shear  layer  (SL),  in  which 
^  =  Re~^/2  ^  (Xj  y  ) ,  and  d</dx  >  0  to  avoid  exponential  growth  as  y  -*■  -<*>  in  SL,  which 
is  governed  by  the  pressure-free  boundary  layer  equations 


with 


u  =  dip /By, 


_3u  jhj;  _3u  _  q  3^u 

3x  3x  By  3y2 


(II— 7 . 3b) 


u (x,°°)  =  1,  u(x,-°°)=0,  <(xgep)  =  0. 


( II— 7 . 3c ,d ,e) 


Separation  occurs  at  x  -  xsep.  The  exponential  dependence  mentioned  arises  from 
the  behavior 


4>  'v  -k(x)  +  A(x,y)exp(ic?  (x)y) ,  y  ->  -® 


( II— 7 . 3f ) 


of  the  SL  solution,  where  the  amplitude  A  is  algebraic  in  y.  As  a  result  the  eddy 
E,  0  <  y  <  S(x),  has  the  development 

V  =  Re  $(x,y)  +.  ..+  A(x,y)exp[Re1^g(x,y)  ]+. . .  (II-7.4a) 


where 


Real (g)  <  0  (II-7 ,4b) 

is  required,  but  g,  A  are  functions  to  be  found,  and  the  condition 


g  ■*  0-  as  y  -*•  S(x)~ 


( I 1—7 .4c) 


A3 


M 

- 

y 


A 
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applies,  to  match  with  (II-7.3f).  Substituting  (II-7.4a)  into  the  Navier-Stokes 
equations  we  find  the  governing  equation 


/-  JL 

dtp 

3  \  -  /; 

!i\2  .  1 

3i\2 

\U  3x  " 

3x 

»)  8  =  ll 

iy/  ( 

3x  ) 

( 1 1—7 . 5a) 


for  the  exponent  function  g,  the  right  hand  side  being  a  viscous  contribution, 
whereas  u,  satisfy  Euler's  inviscid  equations 


3tji  u3u  9^  3u  3p  u3^iji  3i(i  3^j  _  _3p 

3y  ’  3x  3x  3y  3x  *  3x^  3x  3x3y  3y  ’ 


( I I— 7 . 5b) 


2  -  -  _ 

or  V  tp  =  f(tp).  In  addition,  near  the  centerline  a  thin  viscous  extra  region  VE  is 
found  necessary  for  the  exponential  terms  (but  not  for  any  algebraic  terms).  In  VE, 
y  =  Re~l/2  y*  with  y*  of  0(1), 


^  =  Re  1  y*uQ(x)+ - +A*(x,y*)exp[Rel^2gD(x) ]+ .  (II-7.6a) 

gQ(x)  =  g(x,o  +),  Real  (gD)  <  0  and  uq(x)  =  u(x,o)  is  the  reduced  centerline  velocity. 
In  VE  the  Navier-Stokes  equations  reduce  to  the  form 


uo(X)  g;  (x) 


3A*  ,  ,  ,  .  33A* 


(II-7.6t) 


for  A*,  where  PQ(x)  is  a  local  pressure  force.  Hence  the  solution  satisfying  the 
symmetry  condition  A*  *  3^A*/3y*^  =  0  at  y*  =  0  is  nontrivial  in  VE, 


.*  */  .  ,  /  */  \  Po<x>  *  *  _  /  ,^1/2  t  ^ 

A  =  a-j^  (x)  sinh(y  (x)y  )  -  ^-y-yy  »  V  =  (uQgg)  (II-7.6c) 

o 

A 

the  function  a^*(x)  being  fixed  in  principle  by  matching  with  E  as  y  +  ®. 

The  behavior  of  SL,  E,  VE  as  x  -*■  00  downstream  is  as  follows.  In  SL  a  momentum 
integral  of  (II-7.3b-d)  gives,  for  all  x  >  xsep 

00 

k(x)  =  /  u(u-l)dy  +  I  (H-7.7) 
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where  Isep  is  the  integrated  momentum  deficit  at  separation,  from  (II-7.3e),  and  so 
is  negative  usually.  That  suggests  that  as  x  “  the  Chapman  form  (fc)  is  attained 
with  a  surprisingly  high  relative  error  of  order  such  that 


1/2 

i|>  ^  x  fc(n)  +  f^(n)+.  ••  (x  -  x) .  (11-7. 8a) 

From  ( 1 1—  7 . 3b-d) ,  fx  satisfies  fi"+  (fc  f1"  +  f c '  f1')/2  =  0  and  f1'  (t«)  =  0. 

Hence 

(r,)  =  Djf’Oi)  +  Ej^  ( 1 1  — 7 .8b) 


where  are  constants  dependent  upon  starting  conditions  and  an  origin  shift 

in  y.  Consistency  with  (II-7.7)  then  yields  the  asymptote 


k(x)vx^^  +  I 


sen 


as 


(11-7 .8c) 


and  E-^  =  Is  .  The  relative  error  in  Kirchhoff's  parabola  (II-7.1)  is  also  generally 
of  the  same  nigh  order,  x”l/2>  in  that 


S(x)M)x^2+  s 


as 


(IT-7 .9) 


with  sQ  ^  0.  This  can  be  established  either  by  analysis  of  the  farfield  in  Kirchhoff's 
free  streamline  solution  or  by  reference  to  the  exact  solution  for  a  broadside-on 
flat  plate  for  instance:  see  page  499  of  Ref.  18. 

The  small  corrections  in  (II-7.8c),  (II-7.9)  control  the  downstream  flow.  Thus 
in  E,  as  x  ->  00  with  a  =  fixed,  tj)  ^  -  kx^^  o/b  +  c^  a  +  0(1)  from  (IT. 7. 5b), 

(II-7.3a),  with  c^  =  (k  aD^/b  -  Isep/b)  positive.  So  (II-7.5a)  requires  the  underlying 
form 

-1/2  -1 

g  ^  gD(o)  +  x  g^ (<?)  +  .  •  .+x  ig2(o)  + -  ( 1 1—7 .10a) 


for  g,  where 
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gD  "  *(o2-b2)/4b,  gx  -  B^a  , 


bBi2  bc^B^  o  c2  2  k  / 

go  =  -  +  -  a  +  Bocj  +  —  a  lna  -  — —  o'  (II-7.10b) 

L  <  k  ^2  16b 


which  satisfies  (11-7. 4b,  c),  provided  B^  =  -k  sQ/2b.  Hence  the  centerline  value 
gQ(x)  =  g(x,o  +) ,  required  for  VE,  has  the  development 

2 

.  gG(x)  'v-  -  <b/4  +  bB^/<x  +. . .  as  x  -*■  °°.  (II-7.10c) 


•jf 

That  leaves  the  exponent  function  y  defined  in  (II-7.6c)  having  the  behavior 

JL  *1 

y  'V'  Bi  x-i  as  x  ■*  so  that  essentially 


A*  «  sinh(B^y*/x) 


(II-7 .10d) 


there.  Further  terms  can  be  investigated  for  large  positive  x,  including  the 
nature  of  E  when  y  remains  finite  and  the  equations  controlling  the  amplitude  A 
in  (II-7.4a).  Also,  as  a  check  we  notice  that  the  exponential  dependence  in  (II-7.4a) 
becomes  exp[<  Re1/2  (y-b  x1/2)/ 2  x^2]  from  (II-7.10a,  b)  as  y  b  x1/2  in  E,  which 
reproduces  the  Chapman  decay  in  SL;  see  (II-7.3f)  when  (II-7.8c)  holds.  The  principal 
concern,  however,  is  that  VE  expands  linearly  «  Re-^2  x  downstream  (Fig.  11),  from 
(II.7.10d),  due  to  the  factor  Bj_,  i.e.,  to  sQ,  i.e.,  to  the  correction  term  in 
(II-7.9).  Significant  viscous  effects  therefore  spread  out  faster  than  the  x^2, 
Re-^2  xl/2  expansion  rates  of  Kirchhoff's  parabola  and  Chapman's  SL  and  are  bound 
to  affect  the  latter  to  some  extent  on  a  longer  scale  downstream.  That  scale  is 
x  =  0  (Re),  from  comparison  of  the  growth  rates  Re-^^2  x  and  x^-/2. 

1  /  2 

At  least  one  adjustment  of  the  solution  dees  occur  first  on  a  scale  x  =  0  (Rex/  ) 
we  note,  where  X  =  Re-!/2  x  is  finite.  There,  in  summary,  VE  expands  to  0(1)  thick¬ 
ness  in  y  and  its  solution  has  the  development 


vp  =  —  KyRe~^2/b  +  . . .  .+  A^(X,y)exp[-Re^2  <b/4]+....  (II-7.11a) 


implied  by  (II-7.10c).  The  Navier-Stokes  equations  yield  the  backward  heat  equation 
-<  32  A^/3  X  3y  *  b  3^  A^/3y^  for  the  amplitude  3A^/3y.  The  solution  satisfying 
symmetry  at  the  centerline  and  the  join  with  both  E  and  VE  as  X  -*■  0  +  and  with  E 
as  y  -*•  00  is  again  nontrivial, 
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a2 

A,  ■  — —  exp 

1  *1/2 


pi)  [ «-  pf) exp  • 


(II-7 . lib) 


Here  a2»  C2  are  constants.  Therefore  downstream,  as  X  -*■  °°,  remains  0(1)  if  y 
is  restricted  to  0(X^/2)  or  less,  but  A^  increases  exponentially  for  y  of  order  X 
and  the  linear  growth  of  VE  noted  earlier  is  maintained  then.  So  the  length  scaling 
x  a-  Re  is  implied  still. 

It  can  be  established  next  that  Kirchhoff's  parabola  (II-7.1)  with  simple 
entraining  motion  underneath  in  E  becomes  unacceptable  on  the  x  =  0(Re)  scale.  For 
suppose  (11-7.1)  does  continue  to  hold,  with  the  motion  in  E  then  having 

ip  =  ip2(X*,Y2)  + - +  A2(X*,Y2)exp[Re1/2  g2  (X*,Y2)  ]  +  . . .  ,  (II-7.12a) 

dominated  by  the  0(1)  entrainment  necessary  for  SL  when  x  =  ReX*  with  X*  finite. 

Here  y  =  Re-*-/2  Y2,  and  0  <  Y2  <  bX*-*-/2  in  E.  From  the  Navier-Stokes  equations  we 
obtain  the  governing  equations 


/  9 _  ^2  3 _ I  /382\2 

\ 2  a?  -  u?  ml 82  =  \7r2 )  • 

(9  9^2  a  \  .  Dg2  .  9g2  9A?  .  9u?  9g?  . 

U2  '  7x*  nj  I  A2  W2  +  U2  SF  +  3X»  2 


( 1 1—7 . 12b) 


9^2  dg?  9a2  9u2  a82  .  _  A /Sg2  \ 

~  3X*  9Y  9Y  ay  9X*  2  |\9Y2  / 


(IT-7 .12c) 


,  Sg2  9^2  a 
3Y2  ’  3Y22  2 


for  g2,  A2,  omitting  certain  unimportant  contributions,  where  u2  =  3il>2 / 9Y 2 *  So 
if  the  simple  entraining  solution  holds  for  all  X*, 


ip 2  ~  ~  <Y2/b,  U2  =  -x/b,  (II-7.13) 
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continuing  the  trend  (II-7.2),  then  (II-7.12b)  gives  g2  ■  <  (o2-b2)/4b  to  merge  with 
(II-7 .10a,  b),  where  a  =  Y2  X*"l/2,  and  (II-7.12c)  reduces  to 


2Y-. 


3Ao 


3Y2 


+  2  X 


9Ao 

— £  +  A, 
9X"  1 


0. 


( I I— 7 . 14a) 


The  general  solution  for  A2  is 


A2  =  X*  1/2  f 2 (Y2/X*) , 


( II— 7 .14b) 


where  the  function  f2  is  to  be  found.  Matching  with  VE,  however,  as  VE  spreads  out 
into  E  (Fig.  11),  demands  that 


f2(C)  11  sinh(c2?) 


( II— 7 .15a) 


from  (II-7.11b).  In  contrast  the  match  with  SL  requires 


f2(S)  =  80b2  exp(-KW2/2b)/C2, 


(II-7 . 15b) 


*  1/2 

from  (II-7.3f)  where  A  ^  60  x  '  for  large  x,  and  the  constant  d2  represents  an  origin 
shift.  The  two  forms  (II-7.15a,  b)  for  finite  5  cannot  be  reconciled  and  so  the 
description  assumed  breaks  down. 

This  breakdown  on  the  0(Re)  length  scale  may  be  taken  as  evidence  that 
Kirchhoff's  parabola  must  be  adjusted  then,  although  other  interpretations  exist. 
Further,  since  an  ellipse  y  =  Re^/2  b  X*l/^  (1-X*/L)^^2  is  the  only  other  eddy 
shape  consistent  with  the  initial  parabola  and  maintaining  enough  uniform  pressure, 
the  conclusion  of  an  0(Re)  closure  length  LRe  follows  as  in  Refs.  1,  2,  14.  These 
references  appeal  instead  to  an  overall  momentum  balance  (as  in  Section  II-5)  to 
suggest  the  closure  length,  whereas  the  breakdown  above  provides  a  perhaps  firmer 
local  cause  for  the  closure. 

The  elliptical  shape  is  found  to  remove  the  contradiction  of(II-7.15a,  b) .  But 
in  turn  it  emphasizes  another  difficult  feature.  It  points  (at  first  sight  -  see 
Section  II-8  below)  to  the  inviscid  closure  problem 


V24i  =  F  (il>) 


( I I— 7 . 16) 
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within  0(1)  distances  of  x  *  LRe.  Here  the  vorticity  F(\p)  is  nonzero,  because  of 
the  incident,  0(1)  thick.  Chapman  profile  in  SL.  So  (II-7.16)  provokes  the  well- 
known  difficulty  concerning  a  strong  0(1)  thick  jet  being  forced  upstream,  out  of 
the  closure  process,  thus  violating  the  slower  flow  assumption  made  for  the 
majority  of  the  eddy  upstream:  see  Fig.  11.  The  backward  jet  is  unavoidable,  to 
preserve  vorticity  along  the  centerline,  in  the  above  account  unless  viscosity  or 
some  other  subtle  mechanism  comes  into  action.  Similar  difficulties  were  discussed 
in  previous  sections  for  the  widespread  cascade.  Proposals  on  this  wake  closure 
difficulty  for  external  flow  are  given  in  Section  8  below.  Meanwhile,  we  turn  to  the 
relation  between  widespread  cascade  flows  and  external  flow,  given  the  importance 
of  the  0(Re)  length  scale  in  the  latter. 

In  the  relation  between  the  widespread  cascade  and  external  flow 
it  is  noted  first  that  while  the  cascade  spread  H  =  Hy/Hc  as  J-AniiL!-*  as  Part  1« 

the  body  scale  flow  and  its  exponential  corrections  are  effectively  as  in  (II-7.3a)  - 
(II-7.6c).  The  difference  is  that  now  S(uo)  =  Hy  is  finite  (Part  I),  instead  of 
(II-7.1)  holding.  Properties  of  the  exponential  terms  then  are  summarized  in  Ref .  17  and 
Appendix  D  (a  novel  aspect  there  is  that  the  exponent  function  g  becomes  complex, 
inducing  oscillatory  behavior).  As  the  cascade  spread  becomes  wide,  for  a  fixed  body 
size  (H  ->•  0) ,  however,  the  parabola  (II-7.1)  emerges  as  a  first  downstream  asymptote, 
encountered  for  1  <<  x  <<  Hc,  before  the  bound  S  S(°°)  reasserts  itself  as  x  becomes 
0(HC).  Also  Hy  =  0(HC1/2)  because  of  (II-7.1).  The  solution  for  the  eddy  shape  S 
during  this  adjustment  is  found  to  be 


TT 


cos-1  {.-*b2x/«%2} 


( I I— 7 .17) 


2  o 

for  large  x  ^  Hy  v  Hc,  matching  with  (II-7.1)  for  small  x/Hy  but  giving  S  -*■  Hy  for 
large  x/Hy2.  Here  0  <  cos--*-  <  ir/2,  while  the  constant  b  is  0(1)  generally,  equal  to 
0.50  for  the  circular  cylinder.  Let  us  write  Hy  =  b^  H^^2  for  large  Hc,  where  b^  is 
of  0(1)  in  general.  Then  H  =  b^  H,.-^2  -*■  0,  and  so  the  widespread  cascade  first  upsets 
the  Kirchhoff  parabola  (II-7.1)  downstream  only  when  x  is  increased  to  0  (H-2) ,  from 
(II-7.17) . 


Comparison  with  the  x  'v  Re  external  closure  scale  indicated  above  therefore 
re-emphasizes  the  stage  H  ^  Re~l/2  or  Hc  Re  as  in  Part  I.  On  the  other  hand, 
the  viscous  closure  for  the  cascade  wake  always  yields  the  scaled  closure  length 
Xreatt»  i*e**  the  unsealed  closure  length  from  Part  1  is  xreatt  =  Re  Hc  - 

Re  bj^  H~^  Xreatf  So  if  Xreatt  'v  Hn,  say,  for  the  widespread  cascade  (H  ■+  0) 
studied  in  Sections  II-2-5,  then  xreatt  %  bj.^  Re  Hn_^.  Three  subcases  arise, 
therefore,  which  can  be  summarized  as  follows.  First,  if  n  =  4,  xreatt  stays  at 
0  (Re)  when  H  -*•  0,  i.e.,  as  (in  a  more  physical  interpretation)  the  cascade  spacing 
Hc  is  increased  for  fixed  Re;  and  when  H  then  falls  through  the  stage  0(Re~l/2) 
there  is  a  continuous  adjustment  from  ( I 1—7 .17)  to  the  elliptical  shape  above 
as  in  Refs.  1,  14.  Second,  if  n  <  4,  xreatt  increases  beyond  0(Re)  when  H  falls, 
at  least  until  the  stage  H  ^  Re"l/2  when  the  closure  first  affects  matters  and 
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reduces  xreatt  back  to  0(Re).  Third,  if  n  >  4,  xreatt  becomes  shorter  than  0(Re) 
as  H  is  decreased,  and  indeed  xreatt  enters  the  adjustment  zone  of  (II-7.17)  when 
H  falls  to  0(Re”l/ (n-2) )  or  hc  is  increased  to  Re^/ (n-2) # 

So,  whatever  the  value  of  n,  the  above  arguments  imply  that  Kirchhoff's  solution 
for  external  flow  does  emerge  from  the  widespread  cascade  solution. 
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II-8 .  NON-ENTRAINING  SHEAR  LAYERS 


Encouraged  by  the  final  remark  of  Section  II-7,  we  seek  an  alternative  account 
for  wake  closure  consistent  with  Kirchhoff's  (body  scale)  flow  but  avoiding  the 
central  difficulty  of  backward  jets  noted  after  (II  — 7.16)  .  Again,  previous  comments 
suggest  that  viscous  forces  must  play  a  key  role,  if  the  objection  (11-7. 16) ff  is 
to  be  placated,  and  that  algebraic  behavior  and  nonuniqueness  within  the  shear 
layer  can  occur.  We  pursue  these  themes  below. 

Section  II-4  (A2)  anticipated  that  algebraic  behavior  can  be  achieved  in  t tic- 
lower  reaches  of  a  detached  shear  layer  near  the  start  of  a  separated  wake,  while 
Section  II-5  showed  the  same  can  occur  at  the  other  extreme,  just  before  wake 
closure.  This  raises  the  issue  of  whether  such  algebraic  dependence  can  persist 
in-between,  for  virtually  all  the  shear  layer's  length,  in  the  widespread  cascade 
or  the  external  flow  in  particular  although  the  application  is  quite  wider.  In 
whatever  context  the  shear  layer  is  governed  by  (II-2.2b-d)  or  (II-7.3b-d), 
involving  essentially  the  same  pressure-free  viscous  boundary  layer  equations  as 
in  the  wake  closure  of  Section  1 1—  5 .  Algebraic  velocity  decay  for  any  x  does 
turn  out  to  be  possible,  as  (say)  y-*  -<*>  in  (11-7. 3b-d),  in  the  form  'v  -  b^  |  y  | 
u  -bj  y | y !  Y— 1  with  y+1  positive  for  decay.  The  governing  equation  (IT-7. 3b)  then 
requires  the  expansion 


4>  %  -b1 1 y  |  ~f  + - +  E(x) 


however. 


with  b^ 


being  constant  and 


(11-8. la) 


E  =  -(y+1)  (r+ 2)  (x-9.)/y 


(11-8. lb) 


where  l  is  an  arbitrary  constant.  So  the  restriction  -1  ■  r  •  1  applies  to 

(II-8.1a)  as  well  as  y  ^  0  to  suppress  any  constant  of  addition  to  tp .  Otherwise, 
if  y  >  1  is  assumed,  (II— 8.1a)  is  replaced  by 

~  -  _  I 

<J/  'v  -6(x-K.)  |  y  |  +...,  v  .  (1  1—8.  If) 

The  leading  term  in  (II— 8.1c)  gives  an  exact  solution  of  the  boundary  layer 
equations,  incidentally,  while  (II —8 .la)  verifies  that  for  the  widespread  cascade 
the  downstream  form  ( 1 1—4 . 1 5b— f)  in  til  merges  satisfactorily  into  the  shear 
layer  I  on  a  longer  length  scale.  The  two  allowable  forms,  (1 1-8. la)  (with 
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0  <  | y |  <  1)  and  (II-8.1c),  present  an  alternative  to  the  classical  exponential 
decay  (see  (II-7.3f))  usually  taken  by  shear  layers  but,  unlike  the  latter, 
(II-8.1a)  do  not  require  any  direct  entrainment  of  fluid  from  below  the  shear 
layer.  This  is  another  encouraging  point,  for  classical  entrainment  tends  to 
ruin  most  previous  models  (Refs.  1,2,13-15)  of  grossly  separated  motion. 

The  non-entraining  shear  layer  can  also  start  satisfactorily  (Fig.  12),  for 
example  by  taking  over  from  an  incoming,  entraining.  Chapman  form  holding  as 

x-K)  +  (as  required  in  (II-2.2d)  for  example),  in  the  following  manner.  The 

**1/2 

relatively  thin  0(xx/^)  Chapman  form  requires  the  entrainment  ^  =  -k  x  '  ,  say 
at  y  =  0  -.  Beneath  this,  therefore,  a  thicker  subzone  is  induced  where  |y|  is 
0(d^(x))  say,  d^  >>  x^^  unknown  and  the  solution  has  the  displaced  form 
of  a  power  law 


4/  =  -<  (-tc^y  +  d^(x))  *  +. . . .  (11-8. 2a) 

where  T,  k^(>0)  are  constants,  d^(x)  >  0,  and  y  <  0.  Here  the  inviscid  form 

(II-8.2a)  is  consistent  with  the  pressure-free  equations  as  x+  0  +  provided 

di  <<  x  or  (T+l)  / r ( T— 1 )  >  0.  But  entrainment  at  y  =  0  -  requires 

—  r  1/2 

-icd^  =  — kx  ,  so  that 


-1/2T 

d  (x)  =  x  .  (11-8. 2b) 

Hence  (11-8. 2a,  b)  apply  for  T  >  1. ^  Consequently  yet  another  subzone  is  provoked, 
as  viscous  forces  enter  play  where  y  falls  to  -0  (x_l/(r-l))  and  ij;  -*•  0(xr^r-^) 
from  (II-8.2a),  say 


;  r/(r-i) 

ip  =  X 


G(0  +.  . 


with  5  =  y/x 


-1/ ( r-i) 


( I 1—8 . 2c) 


and  £  is  finite  and  negative.  Then  G  satisfies  the  viscous  Falkner-Skan  type 
equation  (I 1-5. 2a)  with  K  =  -  T  <  -1  now,  and 


G  'v  — <  k r  (-£)  r  as  C  0  -  (II-8. 2d) 

to  join  with  (II-8.2a).  The  required  solution  of  (II-5.2a),  (II-8.2d)  has  the 
property  (Ref.  17) 


G  ^  -6  (-C) 


-1 


as  C 


—  00 


(11-8. 2e) 
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consistent  with  (II-5.2a),  while  (II-5.3a)  still  hold  provided  the  middle 
integration  range  there  runs  from  positive  to  £  rather  than  0  to  5, 
leaving  negative  vorticity 


G"  <  0  for  all  £  <  0,  ( II-8 . 2f ) 

in  line  with  (II-8.2d,e),  since  r  >  1.  We  note  that,  since  d^  is  large  and  y 
negative,  ip  is  small  in  (II-8.2a),  as  it  is  in  (II-8.2c),  meaning  that  all  the 
motion  in  the  algebraic  subzones  below  the  Chapman  layer  comes  to  rest  as  x-+0  + 
as  required.  Moreover,  for  small  positive  x  reversed  flow  is  forced  to  occur, 
from  (II-8.2a-e),  so  that  the  algebraic  shear  layer  immediately  becomes  subject 
to  the  influence  of  the  conditions  further  downstream.  This  seems  a  desirable 
feature  in  the  flow.  The  emergent  behavior  (11-8. 2e)  for  small  x  is  the  analogue 
of,  and  continues  into,  the  case  (II-8.1c)  for  finite  positive  x,  provided  the 
origin  i  in  (II-8.1c)  is  zero. 

The  above  indicates  that  the  solution  of  the  characteristic  shear  layer 
problem  posed  in  (II-2.2b-d)  is  not  unique,  then.  The  further  aspects  below 
reinforce  that  conclusion. 

First,  another  successful  start  at  x  =  0  +  is  possible  for  the  non-entraining 
shear  layer  when  the  decay  condition  (II-8.1a)  holds  instead  of  (II-8.1c),  or  when 
i  in  (II-8.1c)  is  nonzero.  In  these  cases  (II-8.2a-f),  beneath  the  incident 
Chapman  form,  are  replaced  by  a  simpler  consistent  description  in  which 

A  A  ^ 

•*  P  [y  -  d^(x)/<2]  as  x  0  +  .  (II-8.3a) 

| 

•  ~  A  A  \ 

Here  ipQ  is  an  arbitrary  profile,  except  that  1 pQ  (y)  'v  -<  (-<2  y)  as  y  -*■  0  -,  with 

constants  <2>  ^  >  0,  and  the  displacement  effect  d£(x)  is  small  such  that 


d2(x) 


1/2X 


( II-8. 3b) 


for  the  entrainment  condition.  Also,  the  behavior 


'P  o(y) 


^  (-y) 


-y 


as  y 


( 1 1-8 . 3c) 


must  hold  with  y  -  1  and  b^  then  identifies  either  with  the  coefficient  b^  in 
(11-8. la)  at  finite  positive  x,  for  y  <  1 ,  or  with  6  1  in  (II-8.1c.)  for  y  =  1. 
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The  restriction  y  -  1  avoids  the  setting  up  of  a  thicker  zone  where  (II-8.2c) 
would  apply  and  would  lead  instead  to  (II-8.2e). 

Thus  both  the  algebraic  decay  forms  (II-8.1a,  c)  can  be  initiated  self- 
consistently  (with  i  zero  or  nonzero  in  (11-8. lc)) ,  branching  away  from  the 
incoming,  exponentially  decaying,  Chapman  form  at  x  =  0  +.  The  algebraic 
start  can  be  from  rest  -  see  (II-8.2a-e)  -  if  (II-8.1c)  holds  with  i  zero,  but 
otherwise  a  nonzero  initial  profile  (II-8.3a)  is  required  (Fig.  12),  below 
the  Chapman  layer.  Both  the  starts  (II-8.2a-e),  (II-8.3a-c),  however,  are 
consistent  with  the  earlier  scaled  inviscid  properties  of  Section  II-4  (a2) 
depending  on  the  value  of  y  there.  Again,  all  the  shear  layer  analysis  is 

A 

subject  to  an  unknown  origin  shift  in  y,  due  to  Prandtl  s  transformation, 
which  adds  arbitrariness  to  the  solution,  while  reversed  flow  in  the  non¬ 
entraining  shear  layer  can  readily  be  present  depending  inter-alia  on  the 
signs  of  the  unknown  constants  b^,  1;  e.g.,  it  occurs  simply  if  x  >  Z  in 
(II-8.1c),  since  u  'v  -6  then. 

Second,  the  non-entraining  shear  layer  can  terminate  in  a  physically 
sensible  fashion  after  a  finite  distance  x  =  X2;  it  cannot  continue  for  all 
positive  x,  we  believe.  The  termination  can  involve  a  singularity  as  x  +  Xj  - 

in  which  the  shear  layer  splits  into  two  main  zones,  1,  3,  each  of  0(1)  thickness 

in  y  and  with  0(1)  velocities  u,  separated  by  a  thickening  viscous  zone  2  of 

width  0(x2-x)-m  in  y  where  the  shear  layer  fluid  is  coming  to  rest  (see  Fig.  12) . 

Here  m  >  0  and  zone  2  has 


\p  'v  const.  +  (x2-x)m+1  G(Z),  Z  =  y/(x2~x)  m.  (II-8.4a) 

So  G  satisfies  (II-5.2a)  again,  but  with  <  =  -(m+l)/m  <  -  1  and  -G* 1 ’  replaces 
G'''.  The  range  of  Z  is  finite  here,  however,  with  minus  the  growth  (II-5.4a) 
being  achieved  as  Z  ZQ  ->  a  similar  growth  occurring  as  (  +  ~Z0  +  ,  all 
subject  again  to  an  origin  shift,  and  the  symmetry  condition  (II-5.2b)  can 
effectively  hold  in-between.  So  (II-5.3a)  now  shows  that 

G"  <  0  for  0  <  Z  <  ZQ,  G”  >  0  for  -  ZQ  <  Z  <  0.  (II-8.4b) 

Centered  around  Z  =  +  ZQ  zone  1  occurs  with  an  arbitrary  shifted  profile 

a  a  ~  —jn 

V  ■*  [y  -  Z0(*2  ~x)  J»  x  "*■  x2“>  ( II-8. 4c) 


consistent  with  the  governing  equations  and  satisfying 
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/m+lx 

-  -  - 

^n(y)  ^  const.  -  C^(-y)  as  y  (II-8.4d) 

(oo)  =  1,  (II-8.4e) 

V11  K  1 


to  match  minus  the  growth  (II-5.4a)  from  zone  2  and  to  comply  with  the  outer 
uniform  stream.  Likewise,  centered  around  £  =  is  zone  3  in  which 


V  +  <P12  [y  +  Zc 


(x2  -  x)  m],  X  ->  x2~, 


(II-8.4f) 


and  the  profile  ifi^2  is  arbitrary  apart  from  the  two  matching  constraints 

in+1 

1^12  (y)  %  const.  +  y  m  as  y  “,  (II-8.4g) 

^12  ^  ^  bl  ^  ®s  y  +  -  ®.  (II-8.4h) 

The  latter  follows  from  (II-8.1a)  for  0  <  |y|  <1  but  from  (II-8.1c)  with  b^ 
denoting  6(£-x2)  for  y=l  and  x2  ^  1. 

Alternatively  the  shear  layer  can  finish  in  a  regular  form  as  x  -*  x2  -  [the 
conditions  (H-8.4e,h)  are  then  satisfied  by  a  single  arbitrary  profile]  or, 
if  x2  =  i,  a  multi-structure  akin  to  that  of  (II-8.2a-e)  can  be  set  up.  The 
precise  form  taken  as  x  ->  x2  -  depends  on  the  subsequent  wake  closure  properties 
and  on  the  flow  beneath  the  shear  layer. 

Underneath  the  shear  layer  a  broader  zone  of 
character  is  controlled  by  whichever  of  (II-8.1a, 
however  the  small  pressure  force  almost  certainly 
yielding  governing  equations  of  the  form 


u  -ii  ,  =  -  p*(x) 

3y  3X  *x  9y 

and  the  matching  condition,  from  (II-8.1a,  c). 


slower  motion  is  produced  whose 
c)  applies.  In  either  case 
reasserts  its  influence, 


+  A 

*1  -2 
9y 


(II-8. 5a) 
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ij/  *v»  -'%2  (^  -  y)  ^  as  y  ■>  6(x)  (II-8.5b) 

Here,  for  (II-8.1a),  X^  =  0  (inviscid  properties,  X2  =  b^  (constant)  and 

0  <  j y |  <  1) ,  whereas  for  (II-8.1c)  X^  =  1  (viscous  properties,  X2  =  6(x-£)  and 

Y  =  1) .  Also,  y  =  3(x)  gives  the  reduced  position  of  the  non-entraining  shear 

layer.  At  the  initial  station  x  =  0  +,  or  at  the  terminal  one  x  -  X£  -,  the 

solution  of  (II-8.5a,  b)  can  start  either  from  an  arbitrary  velocity  profile  u 

(this  is  appropriate  to  Section  II-4  (a2)  for  the  widespread  cascade  for  instance. 

where  u,  y,  y  above  stand  for  0,  Y,  YQ  respectively,  and  the  X-scaling  is 

X  -*■  Y+3) / (Y+l)x)  or  from  a  similarity  form.  For  ^he  latter,  if  X]_  =  1  and 

x  -*  X2~  say,  we  have  ip  %  (x2~x)'*'-M  G(0 ,  ?  =  y/(x2~x)  ,  so  that  (II-5.2a)  holds 

but  with  K=  (1-M) /M  and  -G'"  replacing  G'*'.  Here  I K I  <1  and  the  £-range  is 

_  M  ' 

finite,  0  i  implying  that  6(x)  ^  £, 2  (x2~x)  if,  as  we  assume  here,  the 

symmetry  condition  is  achieved  directly  in  this  broader  zone,  i.e.  (II-5.2b)  also 
applies.  We  would  expect  M  >  0,  for  the  onset  of  closure  (6  -+  0) ,  and  M  >  1/2, 
for  the  scaled  velocity  u  to  increase  then.  Hence  -1  <  K  <  1.  Solutions  of 
(II-5.2a,  b)  then  exist  (Ref.  17)  and  yield  the  required  property 


G  *  -6  U2  -  C)"1  as  K  -  e2- 


(II-8.5c) 


which  satisfies  (II-8.5b).  Solutions  are  also  possible  if  the  pressure  p  remains 
influential  as  x  -►  X2  see  Ref.  17  and  Appendix  E. 

In  summary,  it  seems  clear  that  the  motion  within  and  beneath  the  non-entraining 
shear  layer  can  remain  self-consistent,  from  its  start  to  its  terminal  form.  So 
the  main  difficulty  remaining  is  to  account  self-consistently  for  the  wake 
closure,  near  or  beyond  x  =  X2«  The  correct  account  there  is  not  so  clear  yet. 

2  9 

It  is  tempting  to  appeal  first  to  (II-7.16),  or  its  counterpart  3  ip/dy*  =  F(p) 
for  thin  layers,  and  conclude  that  a  finite  fraction  (where  ip  <  0)  of  the  shear 
layer  should  be  reversed  and  turned  upstream  during  closure  (Fig.  12),  due  to 
local  inviscid  and  pressure  action.  Certainly  the  splitting  in  (II-8.4a-h)  seems 
to  herald  this  and  the  strong  backward  jet  so  produced  would  continue  to  entrain 
no  fluid,  which  is  encouraging.  Each  of  the  conditions  (II-8.1b,  c) ,  however, 
requires  a  linear  decrease  of  the  velocity  contribution  0(y-^)  whose  continuity 
around  the  eddy  therefore  cannot  be  reconciled  readily  with  smooth  inviscid 
turning  near  closure  at  one  end  and  with  the  slower  flow  or  smooth  turning 
necessary  back  towards  the  body  scale  flow  at  the  other  end. 

It  is  tempting,  secondly,  to  imagine  the  lower  part  (II-8.4f-h)  of  the 
terminating  shear  layer  proceeding  smoothly  through  x  =  X2  with  no  significant 
velocity  or  pressure  change  and  coming  out  as  a  starting  profile  for  the  viscous 
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closure  studied  in  Section  II-5.  There,  relatively  near  the  centerline,  the 
requirement  |k|  >  1  leads  to  (II-5.4a),  demanding  in  effect  that  (II-8.1c) 
holds  and  £  =  X£  is  the  closure  length.  That  is  all  consistent  so  far,  but 
back  at  x  =  0  +,  the  initial  shear  layer  then  has  the  form  (II-8.3a-c),  with 
b^  =  6  £,  y  =  1,  and  requires  an  0(1)  velocity  profile  below  the  incident 
Chapman  form.  The  origins  of  such  a  lower  0(1)  profile  are  hard  to  explain. 

A  third  temptation  is  to  postulate  that  a  strong  backward  centerline  jet 
with  edge  decay  different  from  the  terminal  shear  layer's  emerges  from  the 
unified  viscous  closure  of  Section  II-5  as  x^_  -*■  0  +  there.  Yet  the  difference 
in  decay  rates  then  requires  a  similarity  solution,  i.e.  (II-5.2a),  to  hold 
between  the  jet  and  shear  layer  for  small  x-^ .  This  forces  (11-5. 4a)  with 
| K [  >1  which  in  turn  implies  that  X£  =  £  for  both  the  shear  layer  and  the 
jet.  So  again  continuity  of  the  velocity  profile  around  the  circuit  of  the 
eddy  appears  likely. 


Thus  the  x-dependence  in  (II-8.1b,  c)  poses  certain  global  difficulties, 
just  as  entrainment  does  in  classical  shear  layers.  However,  a  novel  description 
of  the  closure  near  and  beyond  x  =  X2  is  possible  now.  Suppose  that  the 
viscous  splitting  (II-8.4a-h)  is  followed  by  the  inviscid  process  of  ( 1 1— 7 . 16) 
locally  again,  but  all  the  streamlines  turned  back  there  remain  part  of  the 
lower  incident  layer  (II-8.4f-h).  This  can  happen  in  principle  because  of 
the  small  velocities  and  reversed  flow  attainable  at  the  edges  of  the  lower  layer, 
in  (II-8.4g,  h) :  it  cannot  happen  with  classical  shear  layers  because  the  latter 
have  at  least  one  nonzero  edge  velocity  with  forward  flow.  (Consider  (I 1-7. 16) 
and  Bernouilli's  law).  So  the  non-entraining  shear  layer  can  "feed  itself",  in 
a  sense,  via  the  closure  process;  see  Fig.  12.  Further  downstream  the  viscous 
closure  of  Section  II-5  can  then  emerge  also.  There  is  no  clear  objection  to  this 


accounts  have  all  failed. 
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II-9 .  FURTHER  COMMENTS 


1.  The  first  clear  sign  to  come  out  of  the  widespread  cascade  analysis  of  Part  II 
is  that  viscous  forces  must  play  a  significant  role  in  massive  wake  closure. 

Inviscid  theory  fails  (Section  II-2) ,  because  of  the  Bernouilli  relation  between 
pressure  and  centerline  velocity. 

That  means  broadly  that  straightforward  estimates  (Refs.  1,2,13-15)  of  eddy 
closure  properties  also  fail  since  they  point  to  predominantly  inviscid  theory. 

Hence  the  explanation  of  massive  wake  closure  must  contain  some  initially  unexpected 
but  realistic  feature(s) . 

The  allied  proposals  of  Section  11-5,8  concerning  non-entraining  shear  layers 
and  viscous  closure  with  algebraic  decay  seem  to  fall  in  line  with  the  above  two 
paragraphs,  as  does  the  long  eddy  result  of  Section  II-3. 

2.  That  non-entraining  shear  layers  can  exist  in  principle,  providing  a  switch 
at  any  stage  from  exponential  entraining  to  algebraic  behavior,  seems  fairly 
definite  (Section  II-8) .  Two  aspects  need  further  examination  nonetheless. 

First,  at  what  stage,  on  what  length  scale,  is  the  switch  made  (if  at  all)? 

For  external  flow.  Section  II-7(a)  suggests  the  x  'v  Re  scale,  in  keeping  with 
most  experimental  and  computational  evidence.  For  the  wide-spread  cascade  there 
is  less  certainty.  Sections  II-3  and  4  (a2)  tend  to  emphasize  lengths  greater 
than  O(H^)  in  X  but  in  retrospect  there  appears  no  obvious  reason  to  really 
discount  shorter  scalings,  despite  the  earlier  estimates.  Indeed,  Section  II-5 
points  to  the  scalings  H4  and  H  (see  near  (II-5 .6a-c) ), Sections  II-2,  4  imply 
that  the  flow  solution  can  no  longer  remain  simple  once  X  reaches  O(H^) ,  while 
the  computations  of  Part  I  seem  to  give  the  H  scaling.  Against  that.  Section  II-3 
shows  that  the  switch  above  may  even  be  unnecessary  for  the  cascade  motion,  in 
which  case  the  non-entraining  shear  layer  can  become  relevant  only  later  in 
external  bluff  body  motion,  as  described  by  Section  II-7,  or  in  other  grossly 
separated  flows. 

This  introduces  the  second  aspect,  namely  that  wake  closure  following  a 
non-entraining  shear  layer  still  seems  to  demand  a  substantial  localized  pressure 
response  (Section  II-8) .  In  external  motion  such  a  response  is  readily  produced 
by  the  blunt  shape  of  the  end  of  the  eddy  (Refs.  1,13  and  Section  1 1  —  7 ) .  By 
contrast  the  widespread  cascade  tends  to  suppress  the  pressure  response  according 
to  both  analysis  (see  Sections  II-2,  4,  5)  and  calculation  (Part  I).  The 
objections  may  well  be  ill-founded,  however,  for  the  alternative  account  of 
Section  II-3  also  induces  a  substantial  pressure,  while  the  upper  portion  of  the 
splitting  in  Section  11-8,  Eqs.  (II-8.4e-e),  could  be  responsible  for  the  required 
pressure  rise  locally. 
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3.  Further  accurate  numerical  study,  of  the  nonlinear  problems  posed  by  (1)  the 
Section  II-3  proposal  for  a  long  eddy  and  branching  solutions,  (2)  the  viscous 
wake  closure  of  Section  II-5,  (3)  the  viscous  reattachment  onto  a  solid  surface 
(Section  II-6),  (4)  the  Section  II-8  non-entraining  shear  layers  of  finite  length, 
and  (5)  the  full  cascade  problem  in  Part  I  for  smaller  values  of  H,  is  desirable 
and  would  add  undoubted  help  to  the  arguments  followed  above  for  massive  stall. 
Otherwise,  such  arguments  always  stay  very  tentative  due  to  the  presence  of 
nonuniqueness,  nonlinearity  and  recirculating  motion.  All  of  ( 1) — ( 5 )  set  not 
inconsiderable  tasks  but  it  is  felt  that  the  prospects  and  insight  into  massive 

stall  properties  provided  by  the  widespread  cascade  problem  are  also  not  inconsiderable 

4.  The  connection  (Section  II-7)  from  the  widespread  cascade  to  external 
bluff-body  motion,  although  hindered  by  our  not  knowing  for  certain  the  power 
governing  the  behavior  Xreatt  «  Hn  as  H-K)  (see  point  2  above),  is  fairly  definite 
on  one  important  feature.  The  Kirchhoff  solution  with  its  parabolic  growth  down¬ 
stream  does  take  control  of  the  body  scale  flow  when  the  cascade  spread  H-^-  is  large. 

This  holds  whatever  the  value  of  n.  The  precise  value  of  n  then  merely  dictates 
what  happens  on  longer  length  scales  downstream  of  Kirchhoff's  parabola.  Given 
that,  it  would  be  very  surprising  if  Kirchhoff's  solution  for  the  body  scale  were 
affected  at  all  significantly  by  further  increase  of  the  cascade  spread  as  external 
flow  properties  take  over  still  more. 

The  sole  hesitation  we  feel  in  stating  the  above  is  slight  and  concerns  the 
possibility  that  the  whole  description  of  wake  closure  for  the  finite-spread 
cascade  in  Part  1  breaks  down  at  a  finite  value  of  H,  allowing  no  limit  form  as 
H-K).  That  would  alter  the  scene  completely.  The  possibility  seems  remote  but 
it  serves  to  stress  the  need  for  more  computational  study  of  the  sensitive  flow 
properties  arising  at  smaller  values  of  H  in  Part  I. 

5.  This  work  has  concentrated  mostly  on  symmetric  wake  closure  with  centerline 
flow.  A  proposal  on  massive  eddy  closure,  or  reattachment,  onto  a  solid  surface 
has  also  emerged  (Section  II-6) ,  however,  which  is  believed  to  merit  further 
attention.  It  suggests  secondary  separation  as  a  possible  consequence  and  finds 
application  in  (e.g.)  the  separating  flow  past  ramps,  past  bodies  with  trailing 
splitter  plates,  and  through  wind  tunnels,  all  of  which  have  attracted  much 
theoretical  and  experimental  interest  in  recent  times.  See  also  Ref.  17.  Again, 
nonsymmetric  wake  closure  is  of  much  practical  and  theoretical  concern. 
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APPENDIX  A 

A  MECHANISM  FOR  UPSTREAM  INFLUENCE 


Consider  an  oncoming  shear  layer  with  virtually  zero  flow  underneath,  which 
seems  not  unreasonable  on  the  verge  of  reattachment.  In  terms  of  an  0(1)  shear  layer 
thickness  for  convenience,  the  oncoming  motion  satisfies 


<p  =  d>Q(y) 
U  =  Uo(y) 


-  -1  cy 

'V'-k  +  a  e  as  y 


ay 

^  e  as  y  +  _<x> 


->  1  as  y  ->■  +  “ 


— oo 


(A-l) 


where  the  constants  k,  a  >  0  due  to  entrainment  upstream.  Then  upstream  influence 
can  occur  on  a  relatively  short  length  scale  A(<<  1)  in  X,  with  a  two-  or  three¬ 
tiered  structure.  In  the  first  tier  where  y  is  0(1),  the  oncoming  form  (A-l)  suffers 
merely  a  displacement  effect,  giving 


*  =  (y  +  A(X) ) ,  U  -  Uq  (y  +  A(X) )  (A-2) 

from  the  controlling  boundary  layer  Eqs.  (1-2. 6a),  where  -A(X)is  the  displacement 
locally,  A(-«)  =  0  and  X  =  AX.  Because  of  the  properties  (A-l)  the  solution  (A-2) 
sets  up  a  second  lower  tier,  at  a  logarithmically  large  distance  below  the  first 
tier,  in  which 

% 

(i|/,U,P(X))  =  l-k  +  A?,  AU,  A2  P(X)]  +  ...  (A-3a) 

and  y  =  a  ^  In  A  +  y,  where  y  is  0(1).  Substituting  again  into  the  boundary  layer 
equations,  we  have  now  the  full  nonlinear  viscous  governing  equations  holding. 


-  dip 

U  =  _  ’ 

0^ 

9U 

dP 

9y 

3X 

9X 

9y 

dX 

with  the  matching  conditions  [from  (A-l),  (A-2)]: 

-  a(y  +  A(X) ) 

U  'V/  e  as  y  +  + 


(A-3b) 


OO 


(A-3c) 
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U  (-»,y)  *  e°y*  P(-“)  =  0 

-  -  -  --  1-2-  (A'3e) 

3i|//3X  (-»,y)  =  -a,  ?(X)  =  -  j  Ufi  (X) 

where  Ug(X)  is  the  slip  velocity  induced  beneath  this  second  tier.  A  third  tier 
may  also  arise,  above  the  first  tier,  to  connect  the  unknown  pressure  P  inviscidly 
with  the  displacement  -A,  but  that  is  secondary  to  the  basic  feature  which  is  the 
behavior  of  the  second  tier  dictated  by  Eqs.  (A-3b-e) .  For  this  admits  upstream 
influence,  as  X  -+  -«,  in  the  form 


iJj  =  -o  X  +  o_1  e°y  +  eGX  f(y)  +  0(e26X) 
P  =  0(e26X) 


(A-4a) 


where  0  is  an  unknown  positive  constant.  From  Eqs.  (A-3b)  with  Eq.  (A-4a) ,  the 
unknown  function  f(y)  satisfies  the  linear  ordinary  differential  equation 


f  ”  '  -of”-  0e°y  (f'  -  a  f)  =  0; 


(A-4b) 


1/2  -1  -  - 

or,  setting  Z  =  2  0  a  exp  (1/2  oy),F=f'  -of,  we  obtain  the  Bessel  equation 


*2 

dZ 


(A-4c) 


for  F(Z).  Hence, 


F(Z)  =  C.  I  (Z)  +  C„  K  (Z)  (A-4d) 

I  o  z  o 

where  C^,  C2  are  constants  and  I0,  KQ  are  Bessel  functions.  However,  IQ  grows  too 
fast  «  Z~l/2  exp  (z)  as  Z  +  1  (y  -*•  »)  to  allow  the  match  (A-3c)  to  be  achieved. 

So  Cj^  *  0.  In  addition  we  require  A  to  be  0  (e^0^)  :  this  is  either  from  the 
linearity  of  the  inviscid  pressure-displacement  law  or  from  the  requirement  of 
negligible  displacement,  A  =  0,  which  often  holds  on  short  length  scales  like  the 
present.  So  |f  ]  <<  eoy  as  y  -*■  +  »,  and  from  Eq.  (A-4d)  we  have  the  solution 

y 

f(y)  =  C0  e°y  f  K  (Z)  e"°y  dy  (A-4e) 
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APPENDIX  B 

FURTHER  UPSTREAM  INFLUENCE 


Appendix  A  described  upstream  influence  arising  in  a  forward-moving  shear  layer 
with  zero  flow  underneath,  but  such  shear  layers  can  also  set  up  a  small  reversed 
motion  underneath  because  of  their  entrainment.  This  Appendix  considers  the  extra 
effects  of  that  small  reversed  motion.  The  main  effect  found  below  is  the  produc¬ 
tion  of  transverse  waves.  These  stem  from  the  local  velocity  and  stream- function 
profiles  of  the  form 


u 


Re'172  -  kl) 

,  _  -1/2  .1  yy  -1 

-k  +  Re  (~  e  -  k^y)  -  Re  yx 


for 


<  y  <  oo 


(B-l) 


with  a  positive  constant,  <  being  the  positive  entrainment  constant  and  y 
being  the  corresponding  positive  index  for  exponential  decay  in  the  shear  layer, 
thus  giving  exponential  growth  in  the  subregion  where  (B-l)  holds,  below  the  shear 
layer.  The  Re  scaling  in  (B-l)  indicates  the  smallness  of  the  velocities  involved, 
although  similar  analysis  in  the  cascade  problem  shows  that  (B-l)  occurs  for  small  X 
instead  there  without  an  Re  scaling. 

The  local  form  (B-l)  is  an  exact  solution  of  the  Navier-Stokes  or  the  boundary 

layer  equations,  interestingly  enough,  related  to  the  asymptotic  suction  solution. 

A  nonlinear  disturbance  to  (B-l)  then  takes  place  if  (ii>,u,p)  =  (-<  +  Re-*^  i, 

-1/2  *  -1  -  1/2  ~  ~  ~ 

Re  u.  Re  p)  within  the  x  =  Re  x  scale  where  y  is  0(1),  so  that  ^,u,p  are 

controlled  by  the  boundary  layer  equations.  Matching  requires  u  'v  exp(yy)  as 

y  -*■  m  and  u  -*■  -k^  as  y  -+  where  k^k^1  =  -p’  for  k-^  nonuniform. 

It  would  be  interesting  to  tackle  the  nonlinear  version  above  since  it  admits 
upstream  influence,  as  does  that  in  Appendix  A.  In  the  present  case  the  upstream 
influence  starts  from  a  small  deviation  *  exp  (Bx)  about  the  exact  solution 
(B-l) .  This  yields  the  linearized  equation 


f’"  -  y  f"  -  (eYy  -  k2)  B  f'  +  y  eYy  6  f  =  a  (B-2) 

for  the  stream  function  perturbation,  with  a  being  a  pressure  constant.  Notice  that 
the  normal  velocity  v  a  -  9tji/3x  plays  a  part  here,  contributing  the  term  in  f V  as 
in  Appendix  A:  (B-2)  is  an  extension  of  (A-4b)  in  fact  incorporating  the  flow 
reversal  due  to  k-^  in  (B-l).  The  solution  of  (B-2)  may  be  found  from  the  substitu¬ 
tions  h  *  Q-^  (f"  -  y  f'),  Q  =  2  y  B^  ^  exp  (yy/2),  which  reduce  (B-2)  to  the 
Bessel  equation 


64 
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for  h(Q)  in  0  <  Q  <  °°,  where  we  have  normalized  k-^  to  be  unity. 

2 

The  solutions  of  (B-3)  are  the  Bessel  functions  I.  (0),  K%I(Q)  where  v  =  1  - 
?  2  >  v  *  v 

4  0  y  •  Assuming  first  that  v  _  0  we  may  discount  the  Iv  function  as  in  Appendix  A 

and  construct  the  form  of  f(y)  by  integration  much  like  that  in  (A-4e) .  The  solu¬ 
tion  decays  exponentially  as  v  4  (Q  ->•  0+)  ,  as  required,  and  so  upstream  influence 
is  admissible,  with  0  -  y  /4.  In  contrast,  if  v  <  0,  so  that  v  =  iv  with  v  real 
and  0  >  y^/4,  then  although  Iv  is  ruled  out  again  by  the  behavior  as  y  4  ®  the 
solution  for  h(Q)  is  oscillatory  a  cos  [v  In  Q]  as  Q  4  0+.  Hence  we  obtain  waves 
a  cos  (vyy/2)  as  y  4  -00.  Upstream  influence  is  still  admissible  with  these  trans¬ 

verse  waves  provided  the  now  significant  motion  set  up  beneath  this  subregion  can 
accommodate  the  waves.  Appendix  C  below  suggests  that  accommodation  is  possible. 
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APPENDIX  C 

SIGNIFICANT  VISCOUS  INFLUENCE  IN  LARGE-SCALE  REVERSED  FLOW 


Throughout  Part  2  the  necessity  of  viscous  action  during  wake  closure  is 
emphasized  (e.g.,  see  Section  2  of  Part  2),  even  though  for  the  wide-spread  cascade 
viscous  forces  seem  largely  negligible  at  first  sight:  again  see  Section  2.  So 
the  question  arises  of  how  viscous  action  may  come  into  play  in  a  predominantly 
inviscid  motion  containing  significant  reversed  flow.  Effectively  this  amounts  to 
considering  the  boundary  layer  equations  for  a  large-scale  flow,  in  the  form 


2 

cM£_  _3_u  9J;  jhj  d£  _1_  9  u 

9y  ’  U  9x  3x  9y  dx  Re  2 

9y 


(C-l) 


where  Re  >>  1  but  x,  y,  u,  p  are  generally  0(1)  rather  than  scaled  in  the 
typical  boundary  layer  fashion. 


At  first  sight  the  viscous  contribution  in  (C-l)  should  be  disregarded,  there¬ 
fore,  but  that  then  points  to  the  backward  jet  inconsistency  referred  to  in  Section 
7  of  Part  2  if  wake  closure  is  to  be  described.  However,  in  view  of  the  short- 
scale  oscillatory  responses  from  viscous  effects  found  in  Appendices  B  and  D,  there 
is  an  alternative  account  in  which  viscosity  matters  more,  as  follows.  Suppose  the 
solution  throughout  depends  on  both  the  general  0(1)  scale  in  y  and  on  the  tradi¬ 
tional  shorter  scale  Y  «  Re1^  y  locally.  Then  a  multiple-scales  treatment  is  called 
for,  with 


9_ 

9y 


1/2 

Re  g'(y) 


(C-2) 


where  g(y)  is  generally  of  0(1) .  The  emphasized  role  of  the  viscous  derivative 
9/9Y  here  allows  a  viscous  perturbation  of  small  amplitude  in  <J>  to  produce  a  sizeable 
response  in  u  and  a  large  response  in  the  vorticity  9u/9y.  Thus  if  a  basic  ^-pro¬ 
file  is  disturbed  linearly  in  the  form 


-1/2 

<J»  =  (y)  +  Re  i|>  (x,y,Y)  + 

D  O 


(C— 3) 


then  (C-2)  gives,  with  u  =  i|>  '  , 

B  B 


9ij; 


1  /? 

u  =  [uR(y)  +  g’(y)  — 2  ]  +  Re“i/Z  — 2  +  . . . 


(C-4) 
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1/2  2 

and  the  disturbance  in  u  is  nonlinear.  Further,  3u/3y  is  dominated  by  Re  (g'(y)) 
32*G/3Y2.  So  if  we  set  ub  +  g'3'l'0/3Y  =  g'2  3ii0/3Y  then  (C-l)  yields  the  full  boundary 
layer  equation  for  ij;0,  in  terms  of  x,Y, 


dip 

_ O 

3Y 


d2ip 

_ o 

8xBY 


_ o 

8x 


3^1^  3"^ti> 

_ o _ o 

3Y2  3Y3 


(C-5) 


with  y-dependence  remaining  passive ,  and  the  range  of  interest  is  -»  <  Y  <  »  because 
of  the  short  scale  involved.  The  pressure  gradient  has  no  effect  here. 

'2  + 

Upstream  where  viscous  effects  decay  tp  -*■  o  and  so  we  require  g  =  _  Ug,  from 
the  y-dependence  in  <J*0 »  and  4>Q  -*•  -  Y.  Hence  the  upstream  profile  for  ij/  (x,y,Y) 
is  simply  a  uniform  stream.  The  stream  is  forward  if  uB  >  o,  since  then  g'  =  Ug1//2, 
but  is  reversed  where  uB  <  o  giving  g'  =  (~Ug)  .  In  consequence,  wherever  reversed 

flow  (ug  <  o)  is  present  upstream  wave-like  viscous  behavior  occurs  transversely  with 


<|>  %  -  Y  +  a  e°X  sin(o1>/2  Y  +  b  )  ,  (C-6) 

O  o  o 

from  (C-5)  as  x  -*•  -°°,  where  a  >  o,  aQ,  bQ  are  constants.  The  waves  span  the  entire 

reversed  flow  zone  ug  <  o.  They  become  a  nonlinear  influence  as  x  becomes  finite 
and  offer  an  alternative  to  the  earlier  inviscid  description.  Again,  wherever  there 
is  forward  flow  upstream,  Ug  >  o,  the  transverse  viscous  response  is  exponential 
with 

ax  -o1/2Y 

*o  ^  Y  +  co  e  e  (C— 7 ) 

if  Y  >  o,  where  cQ  is  a  constant. 

Thus  viscous  effects  can  exert  a  significant  influence  even  in  a  large-scale 
separated  flow,  and  despite  the  largeness  of  Re  in  (C-l)  the  flow  properties  need 
not  be  mainly  inviscid.  A  backward  jet  can  then  be  avoided,  in  particular,  if  wake 
closure  takes  place.  Higher  order  terms  overall  still  require  study  but  the  sus¬ 
tained  waves  in  (C-6)  tie  in  with  those  of  Appendices  B  and  D,  while  the  transition 
between  (C-6)  and  (C—  7)  is  achieved  by  means  of  internal  layers  near  where  uB  =  o. 

A  continued  study  of  the  multiple  scales  effects  here  should  be  interesting. 
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APPENDIX  D 

EXPONENTIAL  BEHAVIOR  IN  REVERSED  MOTIONS 

Here,  for  definiteness,  the  reversed  flow  occurring  far  downstream  on  the  body 
scale  (Part  2  -  Section  7),  i.e.,  near  the  start  of  the  viscous  closure  process 
(Part  2  -  Sections  2-4) ,  is  considered  for  the  cascade  of  finite  spread  H.  Non¬ 
uniqueness  corresponding  to  the  influence  of  conditions  further  downstream  has  been 
noted  already  in  Part  1  for  algebraically  higher  order  terms,  but  the  situation 
is  aggravated  by  the  more  sensitive  exponentially  small  corrections  present  as 
follows . 

The  simplest  starting  form  holding  below  the  Chapman  layer,  for  small  X,  has 
the  reversed  flow  given  by  the  inviscid  solution  \|>  ^  -<  X^^  Y  H  ^  (1-H)  for 

0  <  Y  <  H,  from  (2.12a)  of  Part  2.  Exponential  viscous  dependence  then  enters  the 
reckoning  through  the  decay  in  the  lower  reaches  of  the  Chapman  layer,  which  requires 
the  underlying  expression 


-k  X1/2  Y  - 

<l> - 777~  +  •••  +  A(Y) 


H(l-H) 


1/2 


+  .  .  . 


(D-l) 


for  o  <  Y  <  H.  Here,  the  functions  A(Y) ,  g(Y) .  and  constant  m  >  o  are  unknown  but 
real  (g)  <  c  for  decay  as  X  -*  o+.  The  boundary  layer  equations  than  give  m  =  1/2  and 
require  g  to  satisfy  the  nonlinear  form 


g  +  Y  g' 


'2 
g  » 


(D-2) 


where  we  omit  certain  finite  multiplying  factors  from  the  equation,  for  convenience, 
and  need  g(l)  =  o  to  match  with  the  Chapman  layer.  In  standard  form  the  solution  of 
the  differential  equation  (D-2)  can  be  written  parametrically 


-  Y4>  where  3Y  =  24>  +  4> 


-1/2 


(D-3) 


for  Y  <  1  now.  This  keeps  g  real  and  negative  but  only  for  1  >  Y  >  2  .  As 

Y  -*■  2~ =  YQ  from  above,  an  irregularity  arises  since  there  g(YQ)  =  -2“®/3, 
g'(o)  =  are  finite,  but  g"  «  (Y-Y  )-^^  from  (D-3).  Hence  g  becomes 

w  I 

complex  in  Y  <  Y0.  The  smoothing  out  for  the  3/2  irregularity  in  g  at  Y  =  YQ  -, 
plus  the  removal  of  the  irregularity  found  to  occur  in  A(Y)  there,  can  be  shown 
to  be  achieved  by  Airy  functions  in  a  subzone  of  width  X^'^  near  Y  =  Y0.  These 
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introduce  a  degree  of  arbitrariness  into  the  solution  for  g  in  Y  <  Y0  but  Real  (g) 
remains  negative  as  required.  This  is  seen  perhaps  most  readily  from  the  cubic 
equation  for  g(Y)  ,  obtainable  fortunately  in  this  case,  from  (D-3) , 

TOO/  o 

16  g  -  24  Y  g  +  (9Y  -  12Y)  g  +  (Y  -  1)  =  0.  (D-4) 

For  (D-4)  can  be  used  to  show  that  g  cannot  become  pure  imaginary  in  0  <  Y  <  1,  so 
that  Real  (g)  must  remain  negative.  Also,  at  the  centerline  (Y  -*■  0+)  g  acquires 
one  of  the  complex  roots  of  g^  =  1/16  and  so  again  keeps  Real  (g)  <  0. 

As  a  result,  exponential-oscillatory  behavior  is  present  in  the  reversed  motion 
near  the  start  of  the  wake  closure  process,  the  oscillations  being  of  small  amplitude 
but  very  high  frequency  transversely.  Similar  oscillations  are  found  to  be 
induced  with  the  other  starting  forms  of  Section  2  in  Part  1,  with  the  start  of  the 
long  eddy  in  Section  3,  Part  2  and  with  the  downstream  form  of  Section  7  to  which 
(D-l)-(D-4)  are  equivalent.  Ref.  17  covers  the  application  of  Section  6,  Part  2. 
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APPENDIX  E 

VISCOUS  WAKE  CLOSURE  WITH  PRESSURE  FORCING 


An  analytical  example  occurs  if  the  non-entraining  shear  layer  of  Section  8, 
Part  2  approaches  closure  linearly,  6  ^  ?2(x2  “  x)  as  x  -*■  X2  -  with  the  constant 
?2  positive.  The  eddy  flow  beneath  acquires  a  similarity  form  then,  corresponding 
to  K  =  0  in  Section  5,  Part  2,  so  that  the  eddy  velocity  u  is  0(x2  -  x)  and  in¬ 
creasing.  The  local  similarity  equation  is 


G"t 


(E-l) 


from  (5.2a),  including  a  sign  change,  since  x  is  increasing  here,  and  a  pressure 
gradient  “  irQ.  The  range  of  concern  is  0  -  C  <  ?2’  there  is  no  reversed  center- 
line  jet  (5.2b)  apply,  although  symmetry  conditions  can  also  be  supported  within 
0  <  £  <  C2  if  a  jet  is  present.  We  take  (5.2b)  here.  Integration  of  (E-l)  then 
gives  the  implicit  solution 

f  _ dh _  =  M1' 2  (E-2) 

Jh  (h-hQ)1/2  (h2  +  hQh  +  h/  +  3xro)1/2  '3' 

o 

for  h(5)  =  G'(C),  where  hQ  =  h(o)  is  unknown.  Equation  (E-2)  satisfies  the  symmetry 
condition  (5.2b)  but  also  yields  h  -*■  ®  as  £  -*■  ^  where 


00 


O 


(E-3) 


-2 

and  h  •v  6  (f^-O  as  £  -*■  ^  “•  This  reproduces  the  required  behavior  (8.5c),  on 
integration  for  G(C),  verifying  that  (8.5c)  is  attainable  with  or  without  a  pressure 
gradient  acting.  Further  examples  are  noted  in  Ref.  17. 


Figure  3  Wake,  Edge  and  Centerline  Velocities 
(a)  Nomenclature 


UTRC83-13 


83-3  62-25 


nd  Centerline  Velocities 
Influence  for  H  =  1/3 


ike,  Edge  and  Centerline  Velocities 
Reverse  Flow  Region  for  H  =  1/4 


CLOSURE  DISTANCE 


COMPUTATIONAL  RESULTS 

- MEDIUM  GRID 

— —FINE  GRID 


UTRC83-13 


Fl9Ur*9  *****  Swclng  Scalings 


83-3-  62- 1 1 


UTRC83-13 


Yfo-A131  185  LARGE-SCALE  SEPARATION  AND  HAKE  CLOSURE/REATTACHHENT  -  2/2 

'  THE  CASCADE  PROBLEM <U)  UNITED  TECHNOLOGIES  RESEARCH 

CENTER  EAST  HARTFORD  CT  F  T  SMITH  MAY  83  UTRC/83-13 
UNCLASSIFIED  N00014-81-C-0381  F/G  20/4  NL 


‘j*C#4U<XK* 


